NAVAL 

POSTGRADUATE 

SCHOOL 

MONTEREY,  CALIFORNIA 


THESIS 


INVESTIGATION  OF  PROPAGATION  IN  FOLIAGE 
USING  SIMULATION  TECHNIQUES 

by 

Chung  Wei  Chan 
December  2011 

Thesis  Advisor:  David  C.  Jenn 

Second  Reader:  Ric  Romero 


Approved  for  public  release;  distribution  is  unUmited 


THIS  PAGE  INTENTIONALLY  LEET  BLANK 


REPORT  DOCUMENTATION  PAGE 


Form  Approved  OMB  No.  0704-0188 
Public  reporting  burden  for  this  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instruction, 
searching  existing  data  sources,  gathering  and  maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send 
comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information,  including  suggestions  for  reducing  this  burden,  to 
Washington  headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington,  VA 
22202-4302,  and  to  the  Office  of  Management  and  Budget,  Paperwork  Reduction  Project  (0704-0188)  Washington  DC  20503. 

1.  AGENCY  USE  ONLY  (Leave  blank)  2.  REPORT  DATE  3.  REPORT  TYPE  AND  DATES  COVERED 

December  2011  Master’ s  Thesis 


4.  TITLE  AND  SUBTITLE  Investigation  of  Propagation  in  Foliage  using 

Simulation  Techniques 

5.  FUNDING  NUMBERS 

6.  AUTHOR(S)  Chung  Wei  Chan 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

Naval  Postgraduate  School 

Monterey,  CA  93943-5000 

8.  PERFORMING  ORGANIZATION 
REPORT  NUMBER 

9.  SPONSORING  /MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

N/A 

10.  SPONSORING/MONITORING 
AGENCY  REPORT  NUMBER 

11.  SUPPLEMENTARY  NOTES  The  views  expressed  in  this  thesis  are  those  of  the  author  and  do  not  reflect  the  official  policy 
or  position  of  the  Department  of  Defense  or  the  U.S.  Government.  IRB  Protocol  number  N/A 

12a.  DISTRIBUTION  /  AVAILABILITY  STATEMENT 

Approved  for  public  release;  distribution  is  unlimited 

12b.  DISTRIBUTION  CODE 

A 

13.  ABSTRACT  (maximum  200  words) 

In  a  foliage  environment,  radio  wave  propagation  is  subjected  to  fading  on  both  large-scales  and  small-scales  that 
impair  the  quality  and  reliability  of  data  link  transmission.  This  has  implications  in  many  military  applications.  An 
example  is  the  performance  of  communications  links  and  unmanned  aerial  vehicle  radio  links  when  the  ground  forces 
are  operating  in  foliage  environments. 

The  purpose  of  this  research  is  to  evaluate  some  simple  models  for  propagation  of  radio  waves  in  foliage 
using  an  electromagnetic  field  simulation  application. 

The  three  dimensional  (3D)  electromagnetic  field  simulation  application,  CST  Studio  Suite,  was  used  in  the 
modeling  and  simulation  process.  Specifically,  the  CST  Microwave  Studio  module  was  used  to  model  the  forest 
using  dielectric  blocks.  Various  combinations  of  forest  dimensions,  material  dielectric  parameters  and  antenna 
placements  were  simulated  to  obtain  propagation  models  of  radio  waves  in  foliage  environment. 

The  simulation  models  are  compared  to  three  empirical  models  presented  in  the  literature  for  propagation  in 
foliage  environment.  Using  the  simulation  model,  we  examined  the  coverage  diagram  for  a  transmitter  antenna 
immersed  in  foliage.  The  results  show  that  the  proposed  simulation  models  provide  a  rough  approximation  to 
radiowave  propagation  in  an  actual  rainforest  environment.  Based  on  the  simulated  results,  the  path  loss  in  foliage  is 
affected  by  the  forest’s  electrical  characteristics,  the  height  of  the  transmitter  and  the  height  of  receiver. 

14.  SUBJECT  TERMS  Radio  wave  propagation,  Foliage,  Forest,  EM  wave  simulation  15.  NUMBER  OE 

PAGES 

_ ns _ 

16.  PRICE  CODE 

17.  SECURITY  18.  SECURITY  19.  SECURITY  20.  LIMITATION  OE 

CLASSIEICATION  OE  CLASSIEICATION  OE  THIS  CLASSIEICATION  OE  ABSTRACT 

REPORT  PAGE  ABSTRACT 

Unclassified  Unclassified  Unclassified  UU 

NSN  7540-01-280-5500  Standard  Form  298  (Rev.  2-89) 

Prescribed  by  ANSI  Std.  239-18 


1 


THIS  PAGE  INTENTIONALLY  LEET  BLANK 


11 


Approved  for  public  release;  distribution  is  unlimited 


INVESTIGATION  OF  PROPAGATION  IN  FOLIAGE  USING  SIMULATION 

TECHNIQUES 


Chung  Wei  Chan 

Lieutenant  Commander,  Republie  of  Singapore  Navy 
B.Eng(EE),  National  University  of  Singapore,  1998 


Submitted  in  partial  fulfillment  of  the 
requirements  for  the  degree  of 


MASTER  OF  SCIENCE  IN  ELECTRICAL  ENGINEERING 


from  the 


NAVAL  POSTGRADUATE  SCHOOL 
December  2011 


Author:  Chung  Wei  Chan 


Approved  by:  David  C.  Jenn 

Thesis  Advisor 


Rio  Romero 
Seoond  Reader 


R.  Clark  Robertson 

Chair,  Department  of  Eleotrioal  and  Computer  Engineering 


THIS  PAGE  INTENTIONALLY  LEET  BLANK 


IV 


ABSTRACT 


In  a  foliage  environment,  radio  wave  propagation  is  subjeeted  to  fading  on  both  large- 
scales  and  small-scales  that  impair  the  quality  and  reliability  of  data  link  transmission. 
This  has  implications  in  many  military  applications.  An  example  is  the  performance  of 
communications  links  and  unmanned  aerial  vehicle  radio  links  when  the  ground  forces 
are  operating  in  foliage  environments. 

The  purpose  of  this  research  is  to  evaluate  some  simple  models  for  propagation  of 
radio  waves  in  foliage  using  an  electromagnetic  field  simulation  application. 

The  three  dimensional  (3D)  electromagnetic  field  simulation  application,  CST 
Studio  Suite,  was  used  in  the  modeling  and  simulation  process.  Specifically,  the  CST 
Microwave  Studio  module  was  used  to  model  the  forest  using  dielectric  blocks.  Various 
combinations  of  forest  dimensions,  material  dielectric  parameters  and  antenna  placements 
were  simulated  to  obtain  propagation  models  of  radio  waves  in  foliage  environment. 

The  simulation  models  are  compared  to  three  empirical  models  presented  in  the 
literature  for  propagation  in  foliage  environment.  Using  the  simulation  model,  we 
examined  the  coverage  diagram  for  a  transmitter  antenna  immersed  in  foliage.  The 
results  show  that  the  proposed  simulation  models  provide  a  rough  approximation  to 
radiowave  propagation  in  an  actual  rainforest  environment.  Based  on  the  simulated 
results,  the  path  loss  in  foliage  is  affected  by  the  forest’s  electrical  characteristics,  the 
height  of  the  transmitter  and  the  height  of  receiver. 
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EXECUTIVE  SUMMARY 


The  objective  of  this  thesis  is  to  investigate  various  simulation  models  for 
radiowave  propagation  in  foliage  environments.  There  are  relatively  few  attempts  to 
conduct  simulation  studies  on  radiowave  propagation  in  foliage.  The  possibility  of 
reliable  simulation  techniques  bridges  the  gap  between  theoretical  and  field  studies  for 
propagation  studies  in  a  foliage  environment.  As  an  alternative  to  conducting  field 
measurements,  simulation  provides  a  means  to  validate  performance  of  communications 
and  radar  systems  during  their  design  phase. 

It  has  been  postulated  in  [1]  that  the  forest  can  be  represented  by  dielectric  blocks 
for  a  frequency  up  to  100  MHz.  These  dielectric  forest  blocks  are  characterized  by  the 
quantities  of  relative  permittivity  er  and  conductivity  a.  Based  on  [2],  the  value  of  the 
relative  permittivity  Sr  ranges  between  1.01  and  1.065.  The  value  of  a  ranges  between 
10'"^  S/m  to  10'^  S/m.  This  established  a  basis  for  studying  propagation  through  forest  via 
understanding  the  wave  mechanism  in  the  actual  propagation  process.  Hence,  the 
simulation  models  are  constructed  based  on  a  single  dielectric  slab  to  represent  the  forest. 
A  vertical  dipole  antenna  is  used  to  as  a  transmitter  in  the  simulations. 

Three  theoretical  propagation  mechanisms  in  the  transmitting-receiving  path  in  a 
forest  environment  are  described  in  [1].  The  three  mechanisms  are  the  forest  geometric 
optical  (GO)  components  (consisting  of  direct  and  reflected  rays),  the  sky  wave  and  the 
lateral  wave.  For  distances  larger  than  0.5  km,  the  GO  components  are  largely 
attenuated.  For  short  distances  (i.e.  much  less  than  100  km),  the  sky  wave  component  is 
neglected.  Hence,  lateral  waves  are  the  predominant  propagation  mechanism  in  a  forest 
environment  over  relatively  short  distances.  The  GO  components  and  lateral  wave  are 
depicted  in  Figure  1. 
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Figure  1.  Propagation  mechanisms  in  a  forest  (From  [1]). 

In  this  thesis,  the  simulation  results  are  compared  against  three  empirical  models 
to  determine  the  accuracy  of  the  simulation  models.  The  models  are  Tewari’s  model  [3], 
the  Jansky  and  Bailey  model  [4]  and  the  LITU-R  model  [5].  The  models  were  developed 
using  experimental  data  obtained  from  actual  measurements  carried  out  in  forests.  The 
models  are  applicable  for  the  high  frequency  (HF)  and  very  high  frequency  (VHF)  bands 
and  vertically  (v)  polarized  radiowave  transmission.  The  effect  of  the  forest  on  the 
coverage  diagrams  of  a  transmitter  antenna  immersed  in  foliage  is  also  examined. 
Compared  to  free  space,  the  electric  field  (E-field)  for  propagation  through  foliage  has 
additional  components  of  foliage  attenuation  and  phase  delay  in  foliage.  Due  to  the 
constructive  and  destructive  interference  of  the  ground  reflected  wave  with  the  direct 
wave,  the  typical  maxima  and  nulls  of  the  fields  in  a  coverage  diagram  are  still  present 
when  transmitting  in  a  foliage  environment.  The  reflections  from  the  air-foliage 
interface  can  be  neglected. 

CST  Studio  Suite™,  which  is  an  electromagnetic  (EM)  simulation  application,  is 
used  in  this  thesis.  Specifically,  the  CST  Microwave  Studio  (MWS)  module  is  used  to 
simulate  the  radiowave  propagation  in  the  dielectric  forest  block.  The  dielectric  forest 
block  and  dipole  antenna  are  constructed  in  CST.  Several  parameters,  such  as  s,-,  a  and 
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the  height  of  the  transmitter  antenna  ht,  are  varied  in  the  simulation  to  examine  their 
effects  on  the  path  loss  in  foliage. 

From  the  various  simulations,  it  is  observed  that  the  simulated  results  for  sparse 
foliage  (i.e.,  er  =  1.01  and  a  =  10'^  S/m)  do  not  closely  follow  the  three  empirical  models. 
Generally,  for  denser  foliage  (i.e.,  er  >  1.05  and  a  >  10'^  S/m),  both  Tewari’s  model  and 
the  Jansky  and  Bailey  model  are  better  approximations  to  the  simulation  results  than  the 
LITU-R  empirical  model  for  receiver  situated  at  greater  heights.  Generally,  the  LITU-R 
model  gives  a  better  approximation  to  the  simulation  results  for  both  transmitter  and 
receiver  situated  near  the  ground.  While  the  simulated  results  do  not  exactly  match  those 
given  by  the  three  empirical  models,  the  difference  between  the  simulated  result  and 
empirical  result  is  approximately  10  dB  in  most  cases.  Given  the  lack  of  knowledge  of 
the  electrical  properties  {sr  and  o)  of  the  actual  foliage,  assumed  values  of  the  electrical 
properties  are  used  in  the  simulations.  Hence,  it  is  expected  that  there  will  be  differences 
between  the  simulated  and  empirical  results.  Nevertheless,  given  the  observed 
difference,  the  simulation  models  can  only  be  considered  as  a  rough  approximation  to 
actual  foliage  propagation. 

The  effects  of  the  forest  on  the  coverage  diagram  of  a  transmitting  antenna 
immersed  in  foliage  are  also  examined  using  simulation.  The  coverage  diagram  for  a 
transmitting  antenna  can  be  obtained  from  the  electric  field  plots  of  a  CST  simulation. 
The  height  of  the  transmitter  in  foliage  and  the  forest’s  electrical  properties  are  varied  in 
the  simulation  in  order  to  examine  their  effect  on  the  coverage  diagram.  From  the 
simulation,  it  is  observed  that  an  increase  in  the  height  of  the  transmitter  inside  foliage 
results  in  more  maxima  (or  lobes).  In  addition,  the  angular  difference  between  the 
maxima  and  nulls  becomes  smaller  as  the  antenna  height  is  increased  until  near  the 
canopy  of  the  foliage.  As  the  transmitter  is  sited  higher  up  in  the  foliage,  the  path  length 
of  the  reflected  wave  in  foliage  becomes  longer  and  results  in  a  lower  E-field  strength. 
For  a  transmitter  sited  inside  a  denser  forest  with  higher  Sr  and  a  values,  the  lobes  are 
shortened  as  compared  to  ones  sited  in  a  less  dense  forest.  This  is  due  to  the  increased 
attenuation  in  denser  foliage,  which  is  a  result  of  higher  Sr  and  a  values. 


During  the  conduct  of  the  various  simulation  runs,  we  made  several  observations 
on  the  limitations  of  CSX  with  respect  to  foliage  simulation.  The  simulations  are  run  on  a 
personal  computer  equipped  with  Intel  Xeon  Quad-Core  2.53  GHz  processor  and  4  GB 
RAM  memory.  One  main  limitation  on  the  simulation  with  respect  to  total  simulation 
time  is  the  maximum  frequency  used.  An  increase  in  the  maximum  frequency  results  in 
an  exponential  increase  in  number  of  mesh  cells  generated  by  CSX.  This  leads  to  a 
significant  increase  in  simulation  time.  In  order  for  a  simulation  run  to  be  completed 
within  12  hours,  the  frequency  used  should  not  exceed  75  MHz  for  a  dielectric  block  of 
size  200  m  by  100  m  by  25  m. 

It  has  been  shown  that  it  is  feasible  to  model  foliage  using  EM  simulation 
applications.  Using  the  simulation  technique  proposed  in  this  thesis,  a  radio  frequency 
(RF)  engineer  can  model  an  RE  system  (especially  the  antenna)  that  is  being  developed 
for  operations  in  foliage  using  CSX.  The  performance  of  the  system  can  then  be 
examined  during  the  design  phase.  This  allows  the  designer  to  change  the  design  before 
the  actual  system  is  built  and  is  available  for  actual  field  tests. 
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I.  INTRODUCTION 


A,  MILITARY  OPERATIONS  IN  THE  TROPICS 

Tropical  rainforests  dominate  large  land  masses  in  the  equatorial  region.  An 
example  is  Southeast  Asia.  Thus,  military  forces  of  countries  in  these  regions  are  trained 
to  operate  in  such  forest  environments.  In  the  modem  battlefield,  network  centric 
operations,  as  well  as  the  utility  of  unmanned  vehicles  are  key  enablers  for  mission 
success.  Hence,  the  ability  to  find  reliable  communications  channels  for  radio  networks 
and  unmanned  aerial  vehicle  (UAV)  datalinks  underpins  the  successful  execution  of  a 
mission. 

When  operating  in  tropical  rainforest,  the  foliage  provides  a  natural  camouflage 
for  ground  forces.  In  order  not  to  break  cover  and  potentially  expose  one’s  own  location 
to  the  opposing  force,  the  military  hardware,  including  antennas,  are  kept  hidden  under 
foliage  cover.  Military  forces  thus  have  to  transmit  their  communications  signals  from 
under  the  foliage  cover.  At  the  other  end  of  the  channel,  the  receiving  party  could  be  an 
UAV  operating  in  the  airspace  above  a  forest  or  a  friendly  unit  operating  either  outside  or 
within  the  forest.  A  typical  scenario  of  military  units  operating  in  a  forest  environment  is 
shown  in  Figure  1.  The  implication  of  such  an  operational  scenario  is  that  the 
transmission  path  taken  will  be  through  the  forest  medium.  Likewise,  any  signal  received 
by  forces  operating  under  foliage  has  to  travel  through  the  same  environment.  The  ability 
to  understand  propagation  mechanisms  in  forest  conditions  will  lead  to  the  design  of 
better  systems  for  operations  in  a  forest  environment. 


Figure  1.  An  example  of  battlefield  communications  application  in  forest 

environment  (From  [1]). 
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1,  Physical  Characteristics  of  Tropical  Rainforest 


The  physical  structure  of  a  typical  tropical  rainforest  is  described  in  this  section. 
A  rainforest  is  divided  into  five  different  layers  [2],  A  graphical  depiction  of  the 
rainforest  profile  is  shown  in  Figure  2.  The  emergent  layer  contains  a  small  number  of 
very  large  trees  which  grow  above  the  canopy  layer,  reaching  heights  of  45-55  m.  The 
canopy  is  the  primary  layer  of  the  forest  and  consists  of  the  leafy  top  of  tall  trees,  which 
forms  a  roof  over  the  remaining  layers.  Shorter  trees  and  tall  shrubs  form  the  understory 
layer.  The  plants  in  this  area  seldom  grow  to  3  m.  The  shrub  layer  and  forest  floor  are 
very  dark.  Few  plants  grow  in  this  area  as  a  result.  Hence,  vegetation  under  the  canopy 
layer  is  relatively  sparse. 
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Figure  2.  Profile  of  a  tropical  rainforest  (From  [3]). 
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As  its  name  implies,  the  rainforest  experiences  rainfall  throughout  the  year.  In 
equatorial  regions,  rainfall  may  be  year  round  without  apparent  "wet"  or  "dry"  seasons; 
although,  many  forests  do  have  seasonal  rains.  Even  in  seasonal  forests,  the  period 
between  rains  is  usually  not  long  enough  for  the  leaf  litter  to  dry  out  completely.  Forests 
farther  from  the  equator,  like  those  of  Thailand,  Sri  Lanka,  and  Central  America,  where 
rainy  seasons  are  more  pronounced,  can  only  be  considered  "semi-evergreen"  since  some 
species  of  trees  may  shed  all  of  their  leaves  at  the  beginning  of  the  dry  season  [3].  The 
seasonal  change  in  the  physical  make-up  of  the  forest  impacts  the  electrical  properties  of 
a  forest  and,  consequently,  how  electromagnetic  (EM)  waves  propagate. 

2.  Electrical  Properties  of  a  Forest 

Electrical  properties  of  a  medium  are  specified  by  its  constitutive  parameters,  i.e., 
permeability,  permittivity  and  conductivity  [4].  Permeability  //  and  permittivity  s  are 
expressed  by  the  following  complex  quantities 

A  =  (1) 

^  (2) 

where  jUo  =  ^Ti  x  10'^  H/m  and  So  =  8.85  x  10'^^  F/m. 

The  constant  is  called  the  relative  permeability  and  has  a  value  of  one  for 

lossless,  non-magnetic  materials.  The  constant  is  called  the  relative  dielectric 

constant  or  relative  permittivity  [5].  The  imaginary  terms  in  (1)  and  (2)  give  rise  to 
propagation  losses  due  to  attenuation  and  absorption.  (The  time  convention  used  here  is 

Tamir  proposed  in  [6]  that  the  forest  could  be  represented  by  dielectric  blocks. 
This  established  a  basis  for  studying  propagation  through  forest  via  understanding  the 
wave  mechanism  in  the  actual  propagation  process.  This  was  in  contrast  to  propagation 
studies  that  are  based  on  empirically  or  statistically  derived  models.  However,  one  main 
challenge  in  using  dielectric  models  is  the  need  for  prior  knowledge  of  the  forest’s 
electrical  properties. 
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The  forest’s  eleetrical  properties  are  related  to  characteristics  of  the  forest’s 
physical  structure,  such  as  density  of  canopy,  vegetation  type  and  size  of  trees.  Given  the 
diverse  composition  of  tropical  rainforest,  the  electrical  properties  differ  for  each 
rainforest.  In  some  cases,  the  electrical  properties  for  a  given  forest  changes  according  to 
seasonal  variations  that  affect  the  canopy’s  leaf  density.  Hence,  radiowave  propagation 
in  a  forest  medium  can  be  explained  via  a  simple  propagation  mechanism  that  is 
dependent  on  the  electrical  properties  of  the  forest. 

B,  OBJECTIVES  OF  THESIS 

Current  and  past  research  on  radiowave  propagation  in  foliage  environments  is 
largely  centered  on  either  theoretical  or  field  studies.  While  there  are  numerous 
simulation  studies  on  propagation  in  other  environments  such  as  urban  terrain,  there  are 
relatively  few  attempts  to  conduct  simulation  studies  on  radiowave  propagation  in 
foliage.  This  can  be  attributed  to  the  lack  of  information  on  the  electrical  properties  of 
actual  forests,  which  makes  it  difficult  to  build  an  accurate  model  to  replicate  actual  field 
results  in  a  simulation. 

The  objective  of  this  thesis  is  to  investigate  various  simulation  models  for 
radiowave  propagation  in  foliage  environments.  The  simulation  results  are  compared 
with  empirical  results  to  determine  the  accuracy  of  the  simulation  models.  The 
possibility  of  reliable  simulation  techniques  bridges  the  gap  between  theoretical  and  field 
studies  for  propagation  studies  in  a  foliage  environment.  As  an  alternative  to  conducting 
field  measurements,  simulation  provides  a  means  to  validate  theoretical  models  of 
propagation  and  aid  in  the  design  of  communications  and  radar  systems. 

C.  THESIS  OUTLINE 

There  are  five  chapters  in  this  thesis.  In  Chapter  I,  the  background,  objectives  and 
scope  of  the  thesis  are  presented.  The  theoretical  models  of  propagation  in  foliage  and 
the  various  factors  affecting  performance  of  radiowave  propagation  in  foliage 
environment  are  discussed  in  Chapter  II. 
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In  Chapter  III,  various  EM  field  simulation  applications  which  can  be  used  to 
model  and  simulate  radiowave  propagation  in  foliage  environment  are  presented.  The 
construction  of  the  model  is  also  presented  here.  In  Chapter  IV,  the  simulation 
parameters  and  simulation  results  are  presented.  In  Chapter  V,  a  conclusion  to  the  work 
performed  in  this  thesis  and  recommendations  for  future  work  are  presented. 
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II.  RADIO  WAVE  PROPAGATION  IN  FOLIAGE 


There  are  three  possible  scenarios  involving  a  transmitter-receiver  pair  in  a  forest 
operating  environment.  In  the  first  scenario,  both  the  transmitter  and  receiver  are  sited 
within  the  forest.  In  the  second  scenario,  the  transmitter  is  sited  within  the  forest,  and  the 
receiver  is  sited  outside  the  forest.  The  third  scenario  is  similar  to  the  second  scenario, 
except  that  the  locations  for  the  transmitter  and  receiver  are  switched.  For  scenarios  two 
and  three,  the  propagation  process  can  be  broken  into  two  parts.  One  part  involves 
propagation  in  the  absence  of  vegetation,  while  the  other  part  involves  propagation  inside 
a  forest  environment. 

In  this  chapter,  the  propagation  mechanism  of  radiowaves  in  a  forest  environment 
as  well  as  several  empirical  propagation  models  for  foliage  are  presented  and  discussed. 

A.  RADIOWAVE  PROGATION  MECHANISM  IN  FOLIAGE 

1.  Representing  the  Forest  as  a  Slab 

In  order  to  understand  the  EM  wave  mechanism  in  an  actual  propagation  process, 
it  has  been  proposed  that  the  forest  could  be  viewed  as  a  dissipative  dielectric  slab  [6]. 
The  basic  slab  structure  is  shown  in  Figure  3.  The  dissipative  dielectric  slab  is  assumed 
to  represent  a  forest  with  an  average  tree  height  of  h.  The  transmitter  is  located  at  a 
height  Zo  above  ground.  The  forest  is  characterized  by  the  complex  dielectric  constant 
given  by  [6] 

-  7  —  =  -  760a7l„  (3) 

(0£o 

where  Sr  denotes  the  average  relative  permittivity,  a  denotes  the  averages  conductivity  of 
the  forest  and  Xo  is  the  wavelength  of  a  wave  with  frequency  f.  The  wavelength  Xo  must 
be  sufficiently  large  for  the  representation  of  the  forest  as  a  uniform,  continuous  medium 
to  be  valid  [6].  Based  on  a  reasonably  assumed  average  separation  between  trees  of  1-5 
m  and  that  the  intervening  space  is  usually  filled  with  foliage  and  other  vegetation,  an 


7 


upper  frequency  of  /  equals  to  100  MHz  (or  a  minimum  wavelength  Xq  of  3  m)  seems 
reasonable  [6]. 


z 


Ground 


Figure  3.  Basic  geometry  of  a  forest  dissipative  dielectric  slab  (From  [6]). 

The  forests’  electrical  parameters  Sr  and  a  are  critical  in  describing  the  EM 
propagation  [7].  Unfortunately,  there  is  very  little  data  for  these  two  parameters.  Based 
on  [8],  the  value  of  Sr,  which  was  derived  from  field  measurements,  ranges  between  1.01 
and  1.065.  The  value  for  a  ranges  from  0.01  to  0.160  mmho/m  (equivalent  to  10'^  S/m  to 
1.6  X  10'"^  S/m).  It  was  also  observed  in  [8]  that  a  larger  corresponds  to  a  larger  a.  In 
[6],  the  range  of  Ey  considered  is  from  1.01  to  1.5,  and  the  range  of  a  considered  is  from 
10'  S/m  to  10'  S/m.  A  dense  forest  exhibits  larger  values  for  Ey  and  a,  while  a  thin  forest 
exhibits  smaller  values.  Hence,  it  is  expected  that  the  lower  limits  (i.e.,  Ey  =  1.01  and 
cr  =  10'  S/m)  and  the  upper  limits  (i.e.,  Ey=  1.5  and  cr  =  10'  S/m)  occur  together. 
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From  [4],  for  a  z  polarized  plane  wave  propagating  in  the  x-direction, 


Eix)=2E^e-^^  (4) 

where  Eo  is  a  complex  constant,  and  y  is  the  propagation  constant  and  can  be  expressed  as 

y^a+jP  (5) 

where  a  is  the  attenuation  constant  (Np/m)  and  P  is  the  phase  constant  (rad/m): 


a  =  co' 


lAE 


ll  +  l  —I  -1 


coe 


(6) 


(i-O) 


/dS 


1  + 


+  1 


\0)S  J 


(7) 


The  attenuation  constant  a  affects  the  rate  of  decay  of  the  wave  travelling  in  the  dielectric 
media. 

Three  theoretical  propagation  mechanisms  in  the  transmitting-receiving  path  in  a 
forest  environment  are  described  in  [6].  The  three  components  are  the  forest  geometric 
optical  (GO)  components  (direct  and  reflected  rays),  the  sky  wave  (or  ionospheric  wave) 
component  and  the  lateral  wave  component.  For  distances  larger  than  0.5  km  from  the 
antenna,  the  forest  GO  components  are  largely  attenuated.  The  sky  wave  component  is 
only  significant  for  distances  larger  than  the  one-hop  distance.  In  [6],  a  practical  distance 
of  over  100  km  is  required  for  the  sky  wave  to  be  significant.  Hence,  lateral  waves  are 
the  predominant  propagation  mechanism  in  a  forest  environment  over  relatively  short 
distances.  The  discussion  in  the  following  section  focuses  on  the  GO  and  lateral  wave 
components,  which  are  depicted  in  Figure  4. 
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Figure  4.  Propagation  mechanisms  in  a  forest  (From  [6]). 


2,  Forest  GO  Contribution 


From  [6],  in  the  absence  of  the  ionosphere,  the  first  scenario  reduces  to  that  of  a 
single  half-space  problem  (Figure  4).  The  basic  waves,  which  appear  as  a  direct  ray  and 
a  reflected  ray,  shown  by  the  trajectories  TR  and  TSR,  respectively,  constitute  the  GO 
contribution  for  this  scenario.  The  electric  field  (E-field)  Ep  corresponding  to  these 
waves  may  be  cast  in  the  form  [6] 


= 30I£ 


L 


■  +  /,. 


-JK  \fE''r 


(8) 


where  /  is  the  current  in  the  transmitting  antenna,  £  is  the  length  of  the  transmitting 
antenna,  is  the  distance  travelled  by  the  direct  wave,  is  the  distance  travelled  by  the 
reflected  wave,  and  ko  is  the  wavenumber  of  air  (vacuum)  and  can  be  expressed  as 

K=(o4^o-  (9) 
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From  (8),  it  can  be  seen  that  Ep  is  inversely  proportional  to  and  rr,  i.e.,  the  field 
strength  is  dependent  on  ‘  and  r/' .  The  reflection  from  the  ground  is  neglected. 

3.  Lateral  Wave  Contribution 

If  the  forest  is  a  lossless  medium,  the  additional  diffraction  fields  would  be 
negligible  eompared  to  the  GO  contributions.  However,  given  that  the  forest  is  a 
dissipative  medium,  the  diffraetion  component  has  a  significant  contribution  to  the  field 
[6]. 

The  lateral  wave  is  represented  by  the  trajeetory  TABR  in  Figure  4.  This  wave 
corresponds  to  the  radiation  emitted  at  the  eritieal  angle  of  total  refleetion  6c  given  by 


The  lateral  wave  appears  when  radiation  occurs  from  a  denser  medium  to  a  less  dense 
one  [6].  From  (10),  it  ean  be  seen  that  the  angle  6c  is  defined  for  real  only.  However, 

the  physical  interpretations  are  still  valid  if  the  losses  are  small  (i.e.,  Im(£^^)  «  |), 

and  the  real  part  of  is  implied  in  (10)  [6]. 

When  the  ray  of  the  lateral  wave  from  the  souree  hits  the  air- forest  boundary  (path 
TA  in  Figure  4),  it  is  refraeted  into  the  air  medium,  where  the  ray  travels  tangentially 
along  the  boundary  (path  AB  in  Figure  4)  while  leaking  energy  baek  into  the  forest 
medium  along  the  direetion  of  6c  [6].  Therefore,  some  of  this  leaked  energy  reaehes  the 
observation  point  in  the  forest  (path  BR  in  Figure  4).  From  [6],  it  is  shown  that  the  lateral 
wave  has  the  form 


□  E^  - 1  COS  (f)  cos  Y  +  sin  y) 

(11) 

E^^^  □  E^  [{s^  - 1)  cos  (j)  cos  Y  +  - 1  sin  y] 

(12) 

E^^^  □  Ei^  sin  (f)  cos  y 

(13) 
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where  (j)  is  the  angle  measured  from  the  x-axis,  y  is  the  angle  of  incline  of  the 
transmitting  dipole  antenna  with  respect  to  the  x-axis  and  El  is  given  by 


E 


L 


60/f  e 
s^-1 


A('i+'/v«rhi) 


(14) 


where  ri  is  the  distance  travelled  by  the  lateral  wave  and  s  is  the  total  separation  of  the 
source  and  observation  points  from  the  forest-air  interface  given  by 

s  =  2h-h^-h^-  (15) 


From  (14),  it  can  be  seen  that  the  electric  field  strength  of  the  lateral  wave  is 
dependent  as  compared  to  for  the  GO  waves.  This  is  due  to  the  continuous 
leakage  of  energy  back  into  the  forest  medium  by  the  lateral  wave  across  its  path  AB  [6]. 

Both  the  GO  and  lateral  wave  contributions  to  the  field  observed  in  a  forest 
environment  have  been  explained.  Other  diffraction  components  are  ignored  since  they 
are  of  lower  order  [6]. 

B,  EMPIRICAL  PROPAGATION  MODELS  IN  FOLIAGE  ENVIRONMENT 


Some  empirical  propagation  models  for  foliage  environment  are  explained  in  this 
section.  These  models  were  developed  using  experimental  data  obtained  from  actual 
measurements  carried  out  in  forests.  While  there  are  many  empirical  foliage  propagation 
models  that  have  been  proposed  in  the  past  thirty  years,  the  models  being  discussed  here 
apply  to  the  high  frequency  (HF)  and  very  high  frequency  (VHF)  bands  and  vertically  (v) 
polarized  radiowave  transmission. 


1. 


Tewari’s  Basic  Transmission  Loss  Model  (Tewari) 


In  [9],  an  empirical  propagation  model  based  on  experimental  data  obtained  in  the 
rainforests  of  India  yields  a  predicted  basic  transmission  loss  Lj,  expressed  as 


4  =  -27.57  +  20  log  /  -  20  log 


- +  — ; 


(16) 
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where  /  is  the  frequency  (MHz),  d  is  the  separation  distance  (m),  a  is  the  rate  of 
attenuation  constant  (dB/m),  and  A  and  B  are  constants  evaluated  from  the  measured  data. 

For  vertically  polarized  transmission,  the  values  of  a,  A  and  B  were  determined 
using  the  experimentally  observed  data  given  in  Table  1. 


Table  1 .  Values  of  a,  A  and  B  (From  [9]). 


Frequency  (MHz) 

a  (dB/m) 

A 

B 

50 

- 

0 

1.9170 

200 

0.0125 

0.4989 

1.8358 

500 

0.0135 

0.3658 

0.9040 

800 

0.0140 

0.2661 

0.5331 

A  plot  of  the  basic  transmission  loss  Lb  versus  separation  distance  d  curve  for 
vertically  polarized  transmission  is  generated  using  MATLAB  and  is  shown  in  Figure  5. 
The  MATLAB  code  is  given  in  the  Appendix. 


Figure  5.  Plot  of  basic  transmission  loss  versus  separation  distance  for  vertically 

polarized  transmission. 
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From  Figure  5,  it  can  be  observed  that,  for  a  given  separation  distance,  path  loss 
increases  with  higher  frequency  in  a  forest  environment.  This  is  expected,  since  at  higher 
frequencies,  the  wavelength  is  comparable  to  the  dimensions  of  the  forest’s  constituents 
such  as  leaves  and  branches. 


2,  Jansky  and  Bailey  Model 


From  [10],  the  Jansky  and  Bailey  empirical  model  was  developed  using  measured 
data  obtained  in  the  tropical  forest  of  Thailand.  The  model  is  valid  for  cases  in  which 
both  transmitting  and  receiving  antennas  are  immersed  in  foliage  and  are  separated  by 
distances  between  8  to  1600  meters.  The  model  could  be  used  for  both  horizontal  and 
vertical  polarization.  The  frequency  range  is  between  25  and  400  MHz.  The  loss  Lb  is 
given  by 


4  =36.57 +  20  log/ -20  log 


Ae 


-\6mad 


B 


(17) 


where  /  is  the  frequency  in  MHz,  d  is  the  separation  distance  in  statute  miles,  and  A,  B 
and  a  are  empirical  constants  given  in  Table  2. 


Table  2.  Constants  for  the  Jansky  and  Bailey  empirical  model  (From  [10]). 


Frequency  (MHz) 

Polarization 

a  (dB/m) 

A 

B 

25 

V 

0.0 

0.0 

0.00212 

50 

V 

0.0 

0.0 

0.00106 

100 

V 

0.045 

0.615 

0.000529 

250 

V 

0.050 

0.759 

0.000443 

400 

V 

0.055 

1.02 

0.000523 

Comparing  (17)  with  (16),  we  observe  that  both  empirical  formulas  have  the  same 


form.  However,  their  constants  are  different  since  the  measured  data  on  which  the 
models  are  based  are  taken  from  forests  in  dissimilar  locations. 

A  plot  of  the  basic  transmission  loss  4  versus  separation  distance  d  curve  for 
vertically  polarized  transmission  was  generated  using  MATLAB  and  is  shown  in  Figure 
6.  The  MATLAB  code  is  given  in  the  Appendix. 


14 


Figure  6.  Plot  of  loss  versus  separation  distance  for  Jansky  and  Bailey  empirical 

model. 

From  Figure  6,  it  can  be  observed  that,  for  a  given  separation  distance,  path  loss 
increases  with  higher  frequency  in  a  forest  environment.  This  is  similar  to  that  observed 
in  Tewari’s  model. 

3.  Lateral  ITU-R  (LITU-R)  Model 

In  [1],  an  empirical  model  that  accounts  for  the  excess  foliage  loss  for  lateral 
waves  is  proposed.  The  model  was  modified  from  ITU-R’ s  model  for  VHF  transmission 
in  a  small  forest.  The  VFIF  frequency  band  is  from  30  to  300  MHz.  The  modified  LITU- 
R  model  is  valid  for  VHF  band  transmission  over  a  foliage  depth  of  up  to  5  km.  The 
expression  for  the  LITU-R  model  is 

=  (18) 
where /is  the  frequency  in  MHz  and  d  is  the  separation  distance  in  meters. 
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The  path  loss  is  obtained  by  ineluding  the  effeet  of  foliage  loss  Lutu-r  in  the  plane 
earth  path  loss  model  [1].  Henee,  the  total  pass  loss  in  a  foliage  environment  ean  be 
expressed  as 

L{dB)  =  40  logjo  id)-20  logjo  )  -  20  log  ( 1 9) 

where  d  is  the  separation  distance  in  meters,  ht  is  the  transmitting  antenna  height  in 
meters  and  h,-  is  the  receiving  antenna  height  in  meters.  From  (19),  it  can  be  seen  that 
path  loss  is  also  affected  by  the  antenna  height.  By  increasing  the  transmitting  or 
receiving  antenna  height,  path  loss  is  reduced. 

A  plot  of  the  path  loss  versus  separation  distance  based  on  the  LITU-R  model  is 
generated  using  MATLAB  and  is  shown  in  Figure  7.  The  MATLAB  code  is  given  in  the 
Appendix.  Both  the  transmitting  and  receiving  antenna  heights  are  at  5  m. 


Figure  7.  Plot  of  path  loss  versus  separation  distance  based  on  LITU-R  model. 

From  Figure  7,  it  can  be  observed  that,  for  a  given  separation  distance,  path  loss 
increases  with  higher  frequency. 
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D,  COVERAGE  DIAGRAMS  WITH  TRANSMISSION  THROUGH  FOLIAGE 

In  a  foliage  environment,  it  is  highly  likely  for  both  the  transmitter  and  reeeiver  to 
be  operating  near  the  surfaee  of  the  earth.  In  sueh  a  situation,  multipath  or  multiple 
refleetions  ean  eause  fading  of  the  signal  [4],  A  flat  earth  is  assumed  in  this  discussion. 
The  geometry  of  the  reflections  from  the  ground  in  a  foliage  environment  is  shown  in 
Figure  8.  Reflections  from  the  boundary  between  the  foliage  and  air  are  neglected. 


Figure  8.  Mutipath  effects  from  a  flat  earth  in  foliage. 

In  Figure  8,  the  transmitter  antenna  is  low  gain  and  assumed  to  be  pointing  at  the 
horizon.  The  receiver  is  located  in  free  space  outside  of  the  foliage.  Hence,  the 
transmitter’s  height  ht  is  less  than  the  forest’s  height  h,  while  the  receiver’s  height  hr  is 
greater  than  h.  If  the  transmitter  is  very  near  to  the  ground,  the  angle  if/  is  both  the 
grazing  angle  and  the  elevation  angle.  The  path  lengths  of  the  rays  in  foliage  are  given 
by  lo,  1 1  and  h. 

The  foliage  attenuation  a  and  phase  constant  in  foliage  ji  are  given  by  (6)  and  (7), 
respectively.  The  phase  constant  in  free  space  can  be  expressed  as 

A =0777=-=^-  (20) 

c  ^ 
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The  reference  field  (or  direct  path  in  free  space)  can  be  expressed  as 


E,= 


(21) 


where  G(Qi)  is  the  gain  of  the  transmit  antenna.  Hence,  the  path  gain  factor  (PGF)  F  is 
given  by 


(22) 


where  E(P),  which  is  the  field  at  point  P,  is  expressed  as  (for  the  case  of  point  P  outside 
the  foliage,  i.e.,  hr  >  h) 


E(P) 


Fm 


FL 


-k)^-jPk 


+r 


a/i2 


Rr 


~jPk2  i^\2~h2  ) 


(23) 


where  T  is  the  reflection  coefficient,  Rn  =  R\  +  R2  and  In  =  /i  +  h-  If  the  ground  is 
assumed  to  be  perfect  electrically  conducting  (PEC),  E  has  a  value  of  —1.  For  the  case 
where  point  P  is  inside  the  foliage  (i.e.,  hr  <  h),  E(P)  is  expressed  as 

E(P)  =  ^-iPR,^-aR,  V^(^2)  ^-jPRn^-aR,,  .  ^24) 

R\2 

A  coverage  diagram  can  be  plotted  with  F  shown  as  contours  in  relation  to  d  and 
hr.  An  example  of  a  coverage  diagram  is  shown  in  Figure  9.  The  diagram  is  plotted  for  a 
vertical  Xol2  dipole  located  at  height  20  m  above  ground  in  25  m  high  foliage.  The 
frequency  is  50  MHz,  =  1.01,  cr  =  10  S/m  and  the  reference  distance  is  1000  m. 
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Figure  9.  Example  of  a  eoverage  diagram  for  a  transmitter  loeated  inside  foliage. 


In  Figure  9,  the  contours  represent  the  PGF  for  10  dB,  0  dB,  —10  dB,  —20  dB  and 
—40  dB  relative  to  free  space  at  the  reference  distance.  From  the  diagram,  the  maxima 
and  nulls  of  the  fields  can  be  observed.  This  phenomenon  is  due  to  the  constructive  and 
destructive  interference  of  the  reflected  wave  with  the  direct  wave. 

C.  SUMMARY 

In  this  chapter,  the  mechanism  of  EM  wave  propagation  through  foliage  was 
examined.  Three  empirical  propagation  models  for  foliage  environment  are  also 
discussed  in  this  chapter.  The  effect  of  a  dielectric  forest  block  on  the  coverage  diagram 
of  a  transmitter  in  foliage  was  also  examined.  In  the  next  chapter,  the  simulation  results 
are  analyzed  with  respect  to  the  theory  discussed  in  this  chapter. 
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III.  EM  FIELD  SIMULATION  APPLICATION 


In  this  chapter,  various  EM  field  simulation  applieations  whieh  ean  be  used  to 
model  and  simulate  radiowave  propagation  in  foliage  environment  are  presented.  In 
addition,  the  method  to  eonstruet  the  foliage  model  and  set-up  the  simulation  using  CSX 
Studio  Suite™  is  also  presented. 

A.  REVIEW  OF  SIMULATION  APPLICATION 

There  are  many  EM  simulation  applieations  in  the  market  to  aid  in  the  study  of 
eleetromagnetie  related  problems.  The  CST  Studio  Suite  ™  was  seleeted  for  this 
problem.  The  CST  Studio  Suite™  (CST)  is  developed  by  Computer  Simulation 
Teehnology  AG.  The  eleetromagnetie  simulation  software  CST  STUDIO  SUITE™ 
eomprises  CSX’s  tools  for  the  design  and  optimization  of  devices  operating  over  a  wide 
range  of  frequencies.  CST  consists  of  several  modules.  The  module  that  is  useful  for 
simulating  radiowave  propagation  in  foliage  is  CST  Mierowave  Studio®  (CST  MWS). 
CST  MWS  is  a  speeialist  tool  for  fast  and  aeeurate  three  dimensional  (3D)  EM  simulation 
of  high  frequeney  problems  [10].  CST  MWS  enables  the  fast  and  aeeurate  analysis  of 
high  frequeney  deviees  sueh  as  antennas,  filters,  eouplers,  planar  and  multi-layer 
struetures  [10]. 

Eor  ereation  of  the  computer  aided  design  (CAD)  models,  CST  Studio  Suite™ 
supports  both  options  of  drawing  the  model  within  CST  or  importing  from  an  external 
file.  There  is  a  wide  range  of  primitives  available  to  ereate  eomplex  struetures. 

B,  APPLICATION  OF  CST  STUDIO  TO  SIMULATION  OF  RADIO  WAVE 

PROPAGATION  IN  FOLIAGE 

The  CST  MWS  module  is  used  to  simulate  radiowave  propagation  through 
foliage.  The  transient  solver  in  CST  is  used  for  the  simulation  in  this  thesis.  The 
transient  solver  is  a  general  purpose  3D  simulator.  The  solver  ean  be  used  to  study  the 
fields  propagating  through  a  eomponent.  It  ean  also  deliver  broadband  frequeney  domain 
results  sueh  as  seattering  parameters  (^-parameters)  [10]. 
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1. 


General  CST  MWS  Studio®  Set-up 


a.  General  Simulation  Settings 

In  order  to  optimize  the  simulation  in  CST,  there  are  several  settings  to 
configure  so  that  the  results  are  sufficiently  accurate  for  analysis.  At  the  same  time,  these 
settings  also  ensure  that  the  simulation  does  not  consume  too  many  computing  resources. 

The  mesh  is  set  to  Hexahedral  type.  Under  the  mesh  properties  (Figure 
10),  the  Lines  per  wavelength  is  set  to  10,  which  provides  a  good  compromise  between 
calculation  time  and  achievable  accuracy.  The  Mesh  line  ratio  limit  is  set  to  the  default 
value  of  10. 


Figure  10.  Mesh  properties  settings  in  CST. 

The  boundary  conditions  are  configured  using  the  selections  shown  Figure 

1 1 .  The  Xmin,  Xmax,  Ymin,  Ymax  and  Zmax  boundaries  are  configured  as  ''open  (add 

space/',  which  acts  as  free  space  and  waves  can  pass  through  the  boundary  with  minimal 

reflections.  The  Zmin  boundary  is  configured  as  a  ground  plane  using  "electric  (Et  =  0)". 
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Figure  1 1 .  Properties  of  the  boundary  conditions. 

For  the  units  used  in  the  simulation  (Figure  12),  dimensions  are  defined  in 
meters  and  frequency  is  defined  in  MHz. 


Figure  12.  Configuration  of  units  used  in  CST. 

The  range  of  frequencies  (i.e.,  minimum  frequency  Fmin  to  maximum 
frequency  Fmax)  to  be  used  in  the  simulation  is  defined  using  the  frequency  Range 
Settings  shown  in  Figure  13. 
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Figure  13.  Frequency  Range  Setting  window. 


b.  Definition  of  Variables 

In  CST,  parameters  can  be  defined  for  use  in  the  simulation  model.  The 
parameters  are  defined  in  the  parameter  list  table  on  the  graphic  user  interface  (GUI). 
The  parameters  shown  in  Table  3  are  defined  for  the  model  used  in  CST  to  simulate 
radiowave  propagation  through  foliage. 


Table  3.  List  of  parameters  defined  in  CST. 


Parameter 

Value 

Description 

Units 

1 

H 

varies  according 
to  simulation  run 

Height  of  dielectric  block 

meters 

2 

L 

varies  according 
to  simulation  run 

Length  of  dielectric  block 

meters 

3 

W 

varies  according 
to  simulation  run 

Width  of  dielectric  block 

meters 

4 

c 

Speed  of  light 

meters/second 

5 

f 

varies  according 
to  simulation  run 

Transmission  frequency 

MHz 

6 

lamda 

c/Cf-lO") 

Wavelength 

meters 

7 

a 

varies  according 
to  simulation  run 

Dipole  length 

meters 

8 

r 

0.0005 -lamda 

Dipole  radius 

meters 

9 

xd 

varies  according 
to  simulation  run 

X-coordinate  of  dipole 

meters 

10 

yd 

varies  according 
to  simulation  run 

Y-coordinate  of  dipole 

meters 

11 

zd 

varies  according 
to  simulation  run 

Z-coordinate  of  dipole  and 
represents  the  transmitter  height,  ht 

meters 

12 

d 

0.025 -lamda 

For  mating  discrete  port  to  dipole 

meters 
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2,  Physical  Construction 


Based  on  [6],  the  forest  ean  be  modeled  as  a  single  homogeneous  dielectrie  bloek 
to  represent  the  entire  vegetation  for  frequeney  up  to  100  MHz.  A  half-wave  {XqIT) 
dipole  is  used  as  the  transmitting  antenna. 

a.  One  Layer  Dielectric  Block 

A  single  homogenous  dieleetric  bloek  with  length  (L),  width  (W)  and 
height  (H)  is  construeted  to  represent  the  forest.  The  dieleetric  block  is  constructed  using 
a  Brick  object  in  CST  with  the  properties  shown  in  Figure  14.  The  variables  used  to 
construct  the  dielectric  block  are  obtained  from  the  variable  list  in  Table  3. 


Figure  14.  Properties  of  dielectric  forest  block. 

The  dielectric  block  is  configured  as  a  dielectric  material,  which  has  the 
properties  shown  in  Figure  15  and  16.  The  material  type  is  set  as  Normal.  The  value  for 
Epsilon  (relative  permittivity  Sr)  is  set  according  to  the  test  case  in  each  simulation  run. 
Since  the  forest  is  considered  a  non-magnetic  structure,  the  value  of  Mue  (relative 
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permeability  jiy)  is  set  to  a  value  of  one.  The  electrieal  conduetivity  a  is  set  accordingly 
for  each  simulation  run. 


Figure  15.  Parameters  for  material  property  of  dielectric  block. 
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Figure  16.  Parameters  for  material  property  of  dielectric  block. 

From  the  settings  in  Figures  14,  15  and  16,  the  forest  block  is  constructed 
in  CSX  and  is  shown  in  Figure  17. 
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Figure  17.  Forest  dieleetric  bloek  eonstrueted  in  CSX. 


b.  Half -Wave  Dipole  Antenna 

The  dipole  antenna  is  modeled  as  a  long  thin  rod  with  an  impedanee  of 
73  The  length  of  the  rod  f2  is  set  to  slightly  less  than  half  a  wavelength.  The  rod’s 
length  affeets  the  matehing  of  antenna.  The  length,  whieh  provides  the  optimum 
transmitted  power  (i.e.,  least  refleeted  power),  is  set  during  each  simulation  run.  The 
method  for  obtaining  the  matched  length  is  described  at  the  end  of  this  section.  The 
radius  of  the  rod  is  set  to  0.00052o.  A  discrete  port  is  inserted  at  the  rod’s  center  to  act  as 
the  excitation  source.  In  order  to  construct  the  dipole  antenna,  a  Cylinder  object  with  the 
properties  shown  in  Figure  18  is  used.  A  small  Brick  object  with  properties  shown  in 
Figure  19  is  also  defined.  The  gap  space  for  inserting  the  discrete  port  is  created  by 
subtracting  the  Brick  object  from  the  Cylinder  object  via  a  Boolean  subtract  operation  in 
CSX.  The  variables  defined  in  Table  3  are  used  to  construct  the  dipole  antenna.  The 
material  of  the  dipole  antenna  is  set  as  PEC  with  the  properties  shown  in  Figure  20. 
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Figure  18.  Properties  of  Cylinder  objeet  used  to  eonstruet  dipole  antenna. 


Figure  19.  Properties  of  Brick  objeet  used  to  ereate  space  needed  for  discrete  port. 
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Material  Parameters:  PEC 


Problem  type;  Default 


General 


Thermal  Mechanics  Density 


General  properties 
Material  name: 


PEC 


Type: 


PEC 


Coated  Material... 


Ill 


Color 


Change.. 


r~l  Draw  as  wireframe 

r~l  Draw  reflective  surface 

0  Allow  outline  display 

r~l  Draw  outline  for  transparent  shapes 

I~1  Add  to  material  library 


0% 


T  ransparency  1 0OZ 


! ^ 

Cancel 

Apply 

Help 


Figure  20.  Properties  of  PEC  material  for  dipole  antenna. 


A  discrete  port  is  required  to  feed  the  dipole  antenna  with  power  for  the 
simulation.  In  CSX,  the  radiated  power  from  the  antenna  is  normalized  to  1  W  or  0  dBW. 
A  Discrete  Edge  Port  is  defined  with  the  properties  shown  in  Figure  21 . 
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Figure  2 1 .  Properties  of  Discrete  Edge  Port  used  to  provide  power  to  dipole  antenna. 


From  the  settings  in  Figures  18,  19,  20  and  21,  the  dipole  antenna  of 
length  il  is  construeted  and  shown  in  Figure  22. 


1 

R«frr«nce  73  0h<K 

P<X.V,7)  •.  -Si.  • 


Figure  22.  Dipole  antenna  eonstructed  in  CSX. 
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The  dipole  antenna  has  to  be  matehed  to  ensure  that  most  of  the  energy  is 
transmitted.  When  the  antenna  is  matehed,  the  reflection  coefficient  or  scattering 
parameter  ( 5jj )  is  the  lowest  at  the  frequency  of  interest.  In  order  to  match  the  antenna,  a 
suitable  dipole  length  has  to  be  obtained  for  use  in  the  simulation.  This  is  achieved  by 
using  the  parameter  sweep  functionality  in  CST.  In  the  parameter  sweep  dialog  box,  the 
range  of  values  to  sweep  and  number  of  steps  are  defined.  The  parameter  sweep  function 
produces  a  series  of  5'-parameter  curves  corresponding  to  different  antenna  lengths.  The 
most  appropriate  length  is  chosen  based  on  the  curve  whose  minima  coincides  with  the 
frequency  of  interest.  This  means  that  is  lowest  at  the  transmission  frequency  for  the 
given  length.  Using  Figure  23  as  an  example,  we  see  that  a  dipole  length  12  of  2.52  m 
gives  the  smallest  at  a  frequency  of  50  MHz.  The  selected  length  is  then  used  as  the 
length  of  the  dipole  antenna  for  the  simulation  to  produce  the  EM  fields  in  the  dielectric 
block. 


Figure  23.  results  for  different  dipole  antenna  lengths. 
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3. 


Fields  Monitors 


In  order  to  monitor  the  fields  produeed  by  the  simulation,  E-field  and  Powerflow 
monitors  are  used.  In  CSX,  the  E-field  monitor  stores  the  eleetrie  field  veetors,  while  the 
Powerflow  monitor  stores  the  Poynting  veetor  of  the  EM  field.  For  eaeh  frequeney  to  be 
monitored,  one  E-field  monitor  and  one  Powerflow  monitor  has  to  be  defined.  A  typieal 
property  window  for  a  field  monitor  is  shown  in  Figure  24. 


Figure  24.  Property  window  for  a  field  monitor. 


The  generated  eleetrie  fields  and  Poynting  veetor  ean  be  exported  to  an  ASCII 
file.  When  exporting  the  generated  field  data,  a  resolution  of  step  size  one  is  seleeted  for 
the  X,  y  and  z  eomponents.  This  is  shown  in  Figure  25.  By  exporting  the  data  to  ASCII, 
the  generated  field  data  can  be  manipulated  in  MATFAB  for  analysis  purposes. 
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Figure  25.  Setting  the  resolution  of  field  data  to  be  exported  in  CSX. 

C.  SUMMARY 

In  summary,  the  steps  for  setting  up  CSX  to  perform  the  simulation  for  radiowave 
propagation  in  foliage  were  covered  in  this  chapter.  Xhe  steps  for  constructing  the  forest 
dielectric  block  and  the  dipole  antenna  were  also  covered.  Xhe  results  and  analysis  of  the 
simulation  runs  are  presented  in  the  next  chapter. 
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IV.  SIMULATION  AND  ANALYSIS 


The  various  simulation  runs  performed  in  CSX  to  model  the  radiowave 
propagation  in  foliage  are  presented  in  this  ehapter.  The  analysis  of  the  simulation  runs  is 
also  diseussed  and  presented  in  this  ehapter. 

A,  COMPARISON  WITH  EMPIRICAL  MODELS 

In  this  seetion,  the  objeetive  is  to  validate  the  results  of  the  simulation  model 
against  empirieal  models.  In  this  simulation  model,  the  forest  is  represented  as  a  single 
dieleetrie  bloek.  The  transmitter  is  a  dipole  antenna  that  is  immersed  inside  the  forest. 
The  generated  Poynting  veetor  data  is  then  proeessed  in  MATLAB  to  obtain  a  loss  versus 
distanee  curve.  The  steps  to  construct  the  model  were  described  in  Chapter  III. 

As  the  empirical  models  did  not  specify  the  range  of  values  for  Sr  and  a  of  the 
forest  from  which  their  measured  data  from  obtained,  a  reasonable  assumption  was  made 
on  the  values  to  be  used  in  the  simulations.  Based  on  [6,  8],  the  range  of  er  used  in  the 
simulation  is  from  1.01  to  1.5,  and  the  range  of  a  is  from  10'^  S/m  to  10'^  S/m. 

In  order  to  generate  the  loss-distance  curve,  only  the  field  data  inside  the 
dielectric  block  are  considered.  Since  there  is  no  receiving  antenna  in  the  simulation,  the 
Poynting  vector  at  any  point  in  space  (coordinates  of  x,  y,  z)  represents  the  available 
power  at  a  receiver  located  at  that  point.  Hence,  the  loss  L  can  be  represented  by 

L(dB)  =  1 0  log  (dB)  - 1 0  log  {dB)  (25) 

where  Pt  is  the  transmit  power  in  W  and  Pr  is  the  receive  power  in  W. 

Given  that  the  radiated  power  is  1  W,  Pt  is  0  dBW  and  (25)  can  be  simplified  to 

L{dB)  =  -mogP^{dB)-  (26) 

In  (25),  L  ignores  the  losses  incurred  in  the  receiver  and  represents  the  path  loss  between 
the  transmitter  and  receiver.  The  MATLAB  code  to  plot  the  loss-distance  curves  is 
given  in  the  Appendix.  The  loss-distance  curve  is  used  to  compare  against  various 
empirical  models  to  validate  the  simulation  results  from  CST. 
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As  there  is  no  requirement  to  examine  the  fields  outside  the  dieleetric  bloek  in  this 
section,  the  background  properties  are  configured  as  shown  in  Figure  26. 


Figure  26.  Background  properties  setting. 

Several  variations  of  the  simulation  model  are  executed  to  examine  the  effects  of 
different  variables  on  the  closeness  of  the  simulation  results  to  the  empirical  models.  The 
variables  are  electrical  properties  of  dielectric  block  (i.e.,  s,.  and  a)  and  height  of 
transmitting  antenna.  The  general  observations  for  the  simulation  results  are  discussed  at 
the  end  of  this  section. 
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1. 


Simulation  Results  of  Propagation  in  a  Dielectric  Block 


Test  case  A1  was  constructed  using  the  values  shown  in  Table  4.  The  transmitter 
is  at  a  height  of  5  m  from  the  ground.  In  order  avoid  effects  due  to  reflections  from  the 
forest  edges,  the  transmitter  is  sited  at  a  position  such  that  it  is  at  least  50  m  away  from  all 
the  side  walls  of  the  dielectric  block. 


Table  4.  Values  used  to  set  up  model  for  test  case  A1 . 


Parameter 

Description 

Value 

1 

11 

Length  of  dipole  antenna 

2.54  m 

2 

f 

Frequency  of  interest 

50  MHz 

3 

H 

Height  of  dielectric  block 

20  m 

4 

L 

Length  of  dielectric  block 

200  m 

5 

W 

Width  of  dielectric  block 

100  m 

6 

xd 

x-coordinate  of  dipole  center 

50  m 

7 

yd 

y-coordinate  of  dipole  center 

Om 

8 

zd 

z-coordinate  of  dipole  center 

5  m 

9 

Sr 

Relative  permittivity 

1.065 

10 

a 

Conductivity 

0.000135  S/m 

Based  on  the  values  in  Table  4,  the  simulated  power  flow  plot  of  the  fields  along 


the  x-z  plane  (i.e.,  y  =  0  m)  in  the  dielectric  block  is  shown  in  Figure  27. 


Type 

Monitor 

Component 

Plane  at  y 

Maximum-2D 

Frequency 


Powerflow  (peak) 

Power  (f=50)  [1] 

Abs 

0 

4.67899  dbVA/m''2  at  49.8007  /  0  /  4.925 
50 


Figure  27.  Power  flow  plot  aty  =  0  m  for  test  case  A1 . 
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From  the  exported  data,  the  loss-distanee  eurves  are  plotted  in  MATLAB  for 
different  transmitter  heights  ht  that  are  parallel  to  the  x-y  plane.  The  height  h,-  also 
represents  the  height  of  a  receiver.  The  loss-distance  curves  for  different  h,-  are  shown  in 
Figures  28,  29,  30  and  31.  In  the  figures,  the  Tewari  model  [9],  the  Jansky  and  Bailey 
model  [10]  and  the  LITU-R  model  curves  are  plotted  for  comparison  with  the  simulated 
results  [1]. 

From  Figures  28,  29,  30  and  31,  it  can  be  observed  that  the  loss-distance  curves 
for  the  simulated  results  are  close  to  those  predicted  by  LITU-R  model  (19)  for  low 
receiver  heights  (i.e.,  h,-  =  1-5  m).  As  the  receiver  height  is  increased,  the  loss-distance 
curves  exhibit  a  downward  shift  towards  the  curves  predicted  by  Tewari ’s  model  (16)  and 
the  Jansky  and  Bailey  model  (17). 


Figure  28.  Loss-distance  curves  for  hr  =  1-5  m  (test  case  Al). 
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Figure  29.  Loss-distance  curves  for  h,.  =  6-10  m  (test  case  Al). 


Figure  30.  Loss-distance  curves  for  h,.  =  1 1-15  m  (test  case  Al). 
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Figure  3 1 .  Loss-distance  curves  for  h,-  =  16-20  m  (test  case  Al). 


2,  Effects  of  Different  and  a  values  on  Simulation  Results 

In  test  case  B,  the  effects  of  different  s,.  and  a  values  on  the  simulation  results 
with  respect  to  the  empirical  models  are  examined.  The  set-up  for  the  simulation  is 
similar  to  Table  4.  However,  in  this  test  case,  the  pairs  of  s,.  and  a  values  as  shown  in 
Table  5  are  used  for  each  simulation  run.  In  addition,  different  lengths  are  used  for  the 
dipole  antenna  so  that  the  antennas  are  matched  at  the  transmitting  frequency.  The 
transmitter  antenna  is  sited  at  5  m  above  the  ground. 

Table  5.  List  of  e,-,  a  and  dipole  length  12  values  used  in  test  case  B  simulation  runs. 


Test  Case. 

Sr 

cr  (S/m) 

n  (m) 

B1 

1.01 

0.00001 

2.6 

B2 

1.08 

0.000135 

2.54 

B3 

1.1 

0.000135 

2.54 

B4 

1.5 

0.000135 

2.54 

40 


For  test  case  Bl,  the  simulated  power  flow  plot  of  the  fields  along  the  x-z  plane 
(i.e.,  j  =  0  m)  in  the  dielectric  block  is  shown  in  Figure  32.  The  transmitter  antenna  is 
sited  at  5  m  above  the  ground. 


Type  Powerflow  (peak) 

Monitor  Power  (f=50)  [1] 

Component  Abs 

Plane  at  y  0 

Maximum-2D  4.64894  dbVA/m'^2  at  50  /  0  /  4.925 
Frequency  50 


Figure  32.  Power  flow  plot  atj  =  0  m  for  test  case  Bl. 


The  “test  case  Bl”  loss-distance  curves  for  different  values  of  hr  are  shown  in 
Figures  33,  34,  35  and  36. 


Figure  33.  Loss-distance  curves  for  hr  =  1-5  m  (test  case  Bl). 

41 


In  Figure  33  it  is  observed  that  the  loss-distance  curves  for  the  simulated  results 
approximate  the  LITU-R  model  at  small  separation  distance  from  the  transmitting 
antenna. 


Figure  34.  Loss-distance  curves  for  h,.  =  6-10  m  (test  case  Bl). 

In  Figure  34,  it  is  observed  that  the  loss-distance  curves  for  the  simulated  results 
approximate  the  LITU-R  model  at  distances  of  40-80  m  from  the  transmitting  antenna. 
In  this  test  case,  the  receiver  antennas  are  sited  at  heights  of  6-10  m  from  the  ground. 

In  Figure  35,  it  is  observed  that  the  loss-distance  curves  for  the  simulated  results 
approximate  the  LITU-R  model  at  distances  of  80-100  m  from  the  transmitting  antenna. 
In  this  case,  the  receiver  antennas  are  sited  at  heights  of  1 1-15  m  from  the  ground. 
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Loss,  L  (dB)  Loss,  L  (dB) 


Figure  35.  Loss-distance  curves  for  hr=\  1-15  m  (test  case  Bl). 


-h  =  16  m 


Distance,  d  (m) 


Figure  36.  Loss-distance  curves  for  hr  =  16-20  m  (test  case  Bl). 
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In  Figure  36,  it  is  observed  that  the  loss-distance  curves  for  the  simulated  results 
approximate  the  LITU-R  model  at  distances  of  100-120  m  from  the  transmitting  antenna. 
In  this  case,  the  receiver  antennas  are  sited  at  heights  of  16-20  m  from  the  ground. 

For  test  case  B2,  the  simulated  power  flow  plot  of  the  fields  along  the  x-z  plane 
(i.e.,  j  =  0  m)  in  the  dielectric  block  is  shown  in  Figure  37.  The  transmitter  antenna  is 
sited  at  5  m  above  the  ground. 


Type  Powerflow  (peak) 

Monitor  Power  (f=50)  [1] 

Component  Abs 

Plane  at  y  0 

Maximum-2D  4.61952  dbVA/m'^2  at  50  /  0  /  4.925 
Frequency  50 


Figure  37.  Power  flow  plot  at  j  =  0  m  for  test  case  B2. 


The  “test  case  B2”  loss-distance  curves  for  different  values  of  hr  are  shown  in 
Figures  38  and  39. 

In  Figure  38,  it  is  observed  that  the  loss-distance  curves  for  the  simulated  results 
approximate  the  LITU-R  model  at  separation  distances  up  to  100  m  from  the  transmitting 
antenna.  In  this  case,  the  receiver  antennas  are  sited  at  heights  of  1-5  m  from  the  ground. 
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Loss,  L(dI3)^  ^  ^  _  Loss,  L(dB) 


Figure  38.  Loss-distance  curves  for  hr  =  1-5  m  (test  case  B2). 


Figure  39.  Loss-distance  curves  for  hr  =  16-20  m  (test  case  B2). 
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Comparing  Figure  39  to  Figure  38,  we  observe  that  as  the  height  of  the  receiving 
antenna  increases  (i.e.,  from  1  m  to  20  m),  the  loss  for  a  given  separation  distance 
increases.  At  higher  receiving  antenna  heights  (i.e.,  16-20  m),  the  loss-distance  curves 
are  closer  to  those  predicted  by  Tewari’s  model  and  the  Jansky  and  Bailey  model. 

For  test  case  B3,  the  simulated  power  flow  plot  of  the  fields  along  the  x-z  plane 
(i.e.,  y  =  0  m)  in  the  dielectric  block  is  shown  in  Figure  40.  The  transmitter  antenna  is 
sited  at  5  m  above  the  ground. 
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Type  Powerflow  (peak) 

Monitor  Power  (f=50)  [1] 

Component  Abs 

Plane  at  y  0 

Maximum-2D  4.61328  dbVA/m^2  at  49.7998  /  0  /  4.925 
Frequency  50 


Figure  40.  Power  flow  plot  aty  =  0  m  for  test  case  B3. 


The  “test  case  B3”  loss-distance  curves  for  different  values  of  h,-  are  shown  in 
Figures  41  and  42. 

In  Figure  41  it  is  observed  that  the  loss-distance  curves  for  the  simulated  results 
approximate  the  LITU-R  model  at  separation  distances  up  to  120  m  from  the  transmitting 
antenna.  In  this  case,  the  receiver  antennas  are  sited  at  heights  of  1-5  m  from  the  ground. 
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Figure  41 .  Loss-distance  curves  for  hr  =  1-5  m  (test  case  B3). 


Distance,  d  (m) 


Figure  42.  Loss-distance  curves  for  h,-  =  16-20  m  (test  case  B3). 
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Comparing  Figure  42  to  Figure  41,  we  observe  that  as  the  height  of  the  antennas 
increases  (i.e.,  from  1  m  to  20  m),  the  loss  for  a  given  separation  distance  increases.  At 
higher  antenna  heights  (i.e.,  16-20  m),  the  loss-distance  curves  are  closer  to  those 
predicted  by  Tewari’s  model  and  the  Jansky  and  Bailey  model. 

For  test  case  B4,  the  simulated  power  flow  plot  of  the  fields  along  the  x-z  plane 
(i.e.,  y  =  0  m)  in  the  dielectric  block  is  shown  in  Figure  43.  The  transmitter  antenna  is 
sited  at  5  m  above  the  ground. 


Type  Powerflow  (peak) 

Monitor  Power  (f=50)  [1] 

Component  Abs 

Plane  at  y  0 

Maximum-2D  4.64543  dbVA/m'^2  at  49.8000  /  0  /  4.925 
Frequency  50 


Figure  43.  Power  flow  plot  aty  =  0  m  for  test  case  B4. 


The  “test  case  B4”  loss-distance  curves  for  different  values  of  hr  are  shown  in 
Figures  44  and  45. 

In  Figure  44,  it  is  observed  that  the  loss-distance  curves  for  the  simulated  results 
approximate  the  LITU-R  model  at  separation  distances  up  to  150  m  from  the  transmitting 
antenna.  In  this  case,  the  receiver  is  sited  at  heights  of  1-5  m  from  the  ground. 
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Figure  44.  Loss-distance  curves  for  h,.  =  1-5  m  (test  case  B4). 
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Figure  45.  Loss-distance  curves  for  hr  =  16-20  m  (test  case  B4). 
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Comparing  Figure  45  to  Figure  44,  we  observe  that  as  the  height  of  the  antennas 
inereases  (i.e.,  from  1  m  to  20  m),  the  loss  for  a  given  separation  distanee  inereases.  At 
higher  reeeiving  antenna  heights  (i.e.,  16-20  m),  the  loss-distanee  curves  are  closer  to 
those  predicted  by  Tewari’s  model  and  the  Jansky  and  Bailey  model. 

In  general,  for  denser  forest  (i.e.,  having  higher  Sr  values  of  at  least  1.08  and  a  on 
the  order  10'"^  S/m  and  higher),  the  simulated  results  approximate  those  predicted  by  the 
LITU-R  model  when  the  receiver  is  sited  at  low  heights  (i.e.,  1-5  m)  in  the  simulation. 
At  higher  antenna  heights  (i.e.,  16-20  m),  both  Tewari’s  model  and  the  Jansky  and  Bailey 
model  provide  a  closer  approximation  to  the  simulated  results  than  the  LITU-R  model. 
This  is  expected  as  both  empirical  models  were  proposed  based  on  measured  data 
obtained  in  dense  forest  with  similar  electrical  characteristics  (i.e.,  Sr  values  of  at  least 
1.08  and  a  of  the  order  10'"^  S/m)  as  that  used  in  the  simulation.  Conversely,  it  is 
observed  that  the  simulated  results  based  on  smaller  values  of  Sr  (i.e.,  1.01)  and  a  (i.e.,  10' 
^  S/m)  do  not  closely  match  the  three  empirical  models.  This  is  evidenced  by  comparing 
the  slope  of  the  loss-distance  curves  for  the  simulated  results  to  the  curves  for  the  three 
empirical  models.  The  slope  of  the  curves  for  the  simulated  results  are  much  flatter  as 
compared  to  the  empirical  models.  This  is  because  the  smaller  Sr  and  a  values  represent  a 
less  dense  forest  that  is  very  much  different  from  the  forests  used  to  obtain  both  Tewari’s 
model  and  the  Jansky  and  Bailey  model. 

3,  Effects  of  Different  Transmitter  Heights  on  Simulation  Results 

In  test  case  C,  the  effects  of  different  transmitter  heights  on  the  simulation  results 
with  respect  to  the  empirical  models  are  examined.  However,  in  this  test  case,  different 
transmitter  antenna  heights  ht,  as  shown  in  Table  6,  are  used  in  each  simulation  run.  The 
height,  H,  of  the  dielectric  block  is  30  m.  The  values  of  the  other  variables  needed  to  set¬ 
up  the  simulation  are  kept  constant  for  all  the  simulation  runs  in  this  test  case  and  are 
shown  in  Table  4. 
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Table  6.  List  of  difference  transmitting  antenna  heights  used  in  test  case  C 

simulation  runs. 


Test  Case. 

ht  (m) 

Cl 

2 

C2 

5 

C3 

10 

C4 

20 

C5 

25 

For  test  case  Cl,  the  simulated  power  flow  plot  of  the  fields  along  the  x-z  plane 
(i.e.,  j  =  0  m)  in  the  dielectric  block  is  shown  in  Figure  46.  The  transmitter  antenna  is 
sited  at  2  m  above  the  ground. 


Type  Powerflow  (peak) 

Monitor  Power  (f=50)  [1] 

Component  Abs 

Plane  at  y  0 

Maximum-2D  4.85 1 67  dbVA/m'^2  at  49.8007  /  0  /  1 .925 

Frequency  50 


Figure  46.  Power  flow  plot  at  j  =  0  m  for  test  case  C 1 . 


For  test  case  Cl,  the  loss-distance  curves  for  different  values  of  hr  are  shown  in 
Figures  47,  48  and  49. 

From  Figures  47,  when  the  receiver  is  at  a  height  of  1-5  m,  it  is  observed  that  the 
loss-distance  curves  for  the  simulated  results  are  bunched  together.  The  curves  for  the 
simulated  results  are  closer  to  the  curves  for  the  LITU-R  model  than  the  curves  for  both 
Tewari’s  model  and  the  Jansky  and  Bailey  model. 
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Figure  47.  Loss-distance  curves  for  hr  =  1-5  m  (test  case  Cl). 


Figure  48.  Loss-distance  curves  for  hr  =\  1-15  m  (test  case  Cl). 
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Figure  49.  Loss-distance  curves  for  hr  =  26-30  m  (test  case  Cl). 


From  Figures  48  and  49,  it  is  observed  that,  as  the  receiver  height  increases,  the 
loss-distance  curves  for  the  simulated  results  exhibit  a  downward  shift  and  move  closer 
to  the  curves  predicted  by  Tewari’s  model  and  the  Jansky  and  Bailey  model. 


For  test  case  C2,  the  simulated  power  flow  plot  of  the  fields  along  the  x-z  plane 
(i.e.,  y  =  0  m)  in  the  dielectric  block  is  shown  in  Figure  50.  The  transmitter  antenna  is 
sited  at  5  m  above  ground. 


Type  Powerflow  (peak) 

Monitor  Power  {f=50)  [1] 

Component  Abs 

Plane  at  y  0 

Maximum-2D  4.68 128  dbVA/m^2  at  49.8007  /  0  /  4.925 
Frequency  50 


Figure  50.  Power  flow  plot  at  y  =  0  m  for  test  case  C2. 
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For  test  case  C2,  the  loss-distance  curves  for  different  values  of  hr  are  shown  in 


Figures  51-53. 


Figure  5 1 .  Loss-distance  curves  for  hy  =  1-5  m  (test  case  C2). 

From  Figure  51,  when  the  receiver  is  at  a  height  of  1-5  m,  it  is  observed  that  the 
loss-distance  curves  for  the  simulated  results  are  bunched  together.  The  curves  for  the 
simulated  results  are  closer  to  the  curves  for  the  LITU-R  model  than  the  curves  for  both 
Tewari’s  model  and  the  Jansky  and  Bailey  model. 

From  Figures  52  and  53,  it  is  observed  that,  as  the  receiver  height  increases,  the 
loss-distance  curves  for  the  simulated  results  exhibit  a  downward  shift  and  move  closer 
to  the  curves  predicted  by  Tewari’s  model  and  the  Jansky  and  Bailey  model. 
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Loss,  L  (dB)  Loss,  L  (dB) 


Figure  52.  Loss-distance  curves  for  hr=\  1-15  m  (test  case  C2). 


Figure  53.  Loss-distance  curves  for  hr  =  26-30  m  (test  case  C2). 
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For  test  case  C3,  the  simulated  power  flow  plot  of  the  fields  along  the  x-z  plane 
(i.e.,  j  =  0  m)  in  the  dielectric  block  is  shown  in  Figure  54.  The  transmitter  antenna  is 
sited  at  10  m  above  ground. 


Type  Powerflow  (peak) 

Monitor  Power  (f=50)  [1] 

Component  Abs 

Plane  at  y  0 

Maximum-2D  4.65549  dbVA/m^2  at  49.8007  /  0  /  9.925 
Frequency  50 


Figure  54.  Power  flow  plot  at  j  =  0  m  for  test  case  C3. 


For  test  case  C3,  the  loss-distance  curves  for  different  values  of  hr  are  shown  in 
Figures  55-57. 


Figure  55.  Loss-distance  curves  for  hr  =  1-5  m  (test  case  C3). 
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From  Figure  55,  when  the  receiver  is  at  a  height  of  1-5  m  and  separation  distance 
of  40-110  m  from  the  transmitter,  it  is  observed  that  the  loss-distance  curves  for  the 
simulated  results  are  closer  to  the  curves  for  the  LITU-R  model  than  the  curves  for  both 
Tewari’s  model  and  the  Jansky  and  Bailey  model. 


Figure  56.  Loss-distance  curves  for  hr=\  1-15  m  (test  case  C3). 

From  Figures  56  and  57,  it  is  observed  that,  as  the  receiver  height  increases  from 
1 1-30  m,  the  loss-distance  curves  for  the  simulated  results  are  fairly  similar.  The  curves 
are  closer  to  the  curves  predicted  by  Tewari’s  model  and  the  Jansky  and  Bailey  model 
than  the  LlTU-R  model. 

For  test  case  C4,  the  simulated  power  flow  plot  of  the  fields  along  the  x-z  plane 
(i.e.,  y  =  0  m)  in  the  dielectric  block  is  shown  in  Figure  58.  The  transmitter  antenna  is 
sited  at  20  m  above  ground. 
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Figure  57.  Loss-distance  curves  for  hr  =  26-30  m  (test  case  C3). 


Type 

Monitor 

Component 

Plane  at  y 

Maximum-2D 

Frequency 


Powerflow  (peak) 

Power  (f=50)  [1] 

Abs 

0 

4.65807  dbVA/m^2  at  49.8007  /  0  /  19.925 
50 


Figure  58.  Power  flow  plot  at  j  =  0  m  for  test  case  C4. 
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For  test  case  C4,  the  loss-distance  curves  for  different  values  of  h,.  are  shown  in 


Figures  59-61. 


Figure  59.  Loss-distance  curves  for  hr  =  1-5  m  (test  case  C4). 

From  Figure  59,  when  the  receiver  is  at  a  height  of  1-5  m,  it  is  observed  that  the 
loss-distance  curves  for  the  simulated  results  are  between  the  curves  for  the  LITU-R 
model  and  both  curves  for  Tewari’s  model  and  the  Jansky  and  Bailey  model. 

From  Figures  60  and  61,  it  is  observed  that,  as  the  receiver  height  increases  from 
1 1-30  m,  the  loss-distance  curves  for  the  simulated  results  are  fairly  similar.  The  curves 
are  closer  to  the  curves  predicted  by  Tewari’s  model  and  the  Jansky  and  Bailey  model 
than  the  LITU-R  model. 


59 


Loss,  L(dB)  Loss,L(dB) 


Distance,  d  (m) 


Figure  60.  Loss-distance  curves  for  hr=\  1-15  m  (test  case  C4). 


Distance,  d  (m) 


Figure  61 .  Loss-distance  curves  for  hr  =  26-30  m  (test  case  C4). 
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For  test  case  C5,  the  simulated  power  flow  plot  of  the  fields  along  the  x-z  plane 
(i.e.,  j  =  0  m)  in  the  dielectric  block  is  shown  in  Figure  62.  The  transmitter  antenna  is 
sited  at  25  m  above  ground. 


Component  Abs 

Plane  at  y  0 

Maximum-2D  4.65509  dbVA/m^2  at  49.8007  /  0  /  24.925 
Frequency  50 

Figure  62.  Power  flow  plot  atj  =  0  m  for  test  case  C5. 


For  test  case  C5,  the  loss-distance  curves  for  different  values  of  hr  are  shown  in 
Figures  63-65. 


Figure  63.  Loss-distance  curves  for  hr  =  1-5  m  (test  case  C5). 
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From  Figure  63,  it  is  observed  that,  as  the  receiver  height  increases  from  1-5  m, 
the  loss-distance  curves  for  the  simulated  results  exhibit  a  downward  shift  and  move 
closer  to  the  curves  predicted  by  curves  Tewari’s  model  and  the  Jansky  and  Bailey 
model. 


Figure  64.  Loss-distance  curves  for  hr=\  1-15  m  (test  case  C5). 

From  Figures  64  and  65,  it  is  observed  that,  as  the  receiver  height  increases  from 
1 1-30  m,  the  loss-distance  curves  for  the  simulated  results  are  fairly  similar.  The  curves 
are  closer  to  the  curves  predicted  by  Tewari’s  model  and  the  Jansky  and  Bailey  model 
than  the  LITU-R  model. 

At  any  given  point  inside  the  foliage,  the  field  is  the  sum  of  a  direct  wave  and 
reflected  wave  and  is  given  by  (24).  When  the  receiver  is  in  foliage,  the  point  P  (as 
shown  in  Figure  8)  is  in  foliage,  and  the  path  length  for  the  direct  wave  Rq  is  given  by 

(21) 
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where  Ah  the  difference  between  the  transmitter  height,  ht  (defined  by  the  variable  zd 
in  the  simulation  setup)  and  receiver  height,  hr  and  is  expressed  as 


^h  =  \\-h\■  (28). 

The  path  length  of  the  reflected  wave  Rn  is  given 

R,,=R,+R,-  (29) 


Figure  65.  Loss-distance  curves  for  hr  =  26-30  m  (test  case  C5). 

From  (24),  the  total  field  at  P,  E(P)  is  related  to  Rq  and  hi  addition,  a  higher 
E(P)  results  in  a  lower  path  loss,  L.  The  effect  of  ht  and  hr  on  E(P)  and  L  is  shown  in 
Table  7. 

From  Table  7,  when  the  transmitter  and  receiver  at  sited  near  the  ground,  the  path 
loss  is  the  lowest.  This  is  substantiated  by  Figure  47  that  represents  the  simulation  run 
with  both  transmitter  and  receiver  near  the  ground,  and  the  loss-distance  curves  for  the 
simulation  results  show  the  lowest  L  as  compared  to  other  simulation  runs.  This  occurs 
because  the  direct  and  reflected  waves  are  nearly  in  phase  for  a  vertical  dipole. 
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Table  7.  Effect  of  and  hr  on  E(P)  and  L. 


Case 

ht  and  hr 

Ro 

Ri,  Ri 

Rj2 

E(P) 

L 

1 

High  ht 

High  hr 

Smaller 

Farger  R  j 

Farger  R2 

Farger 

Fower 

Higher 

2 

Fow  ht 

Fow  hr 

Smaller 

Smaller 

Ri 

Smaller 

R2 

Smaller 

Highest 

Fewest 

3 

High  ht 

Fow  hr 

Farger 

Farger  R  j 

Smaller 

R2 

Farger 

Fower 

Higher 

4 

Fow  ht 

High  hr 

Farger 

Smaller 

Ri  Farger 

R2 

Farger 

Fower 

Higher 

For  the  scenarios  described  by  cases  3  and  4  in  Table  7,  the  path  loss  is 
approximately  the  same  since  the  effect  of  a  smaller  Ri  is  negated  by  a  larger  R2  or  vice 
versa.  Similarly  for  case  1,  the  effect  of  a  larger  path  length  for  the  reflected  wave  is 
compensated  by  a  smaller  path  length  for  the  direct  wave  and  results  in  a  path  loss 
comparable  to  cases  3  and  4.  This  is  evidenced  from  the  simulation  results  in  test  cases 
C2,  C3,  C4andC5. 

4,  Observations 

From  the  various  simulation  runs,  it  is  observed  that  the  simulated  results  for 

sparse  foliage  (i.e.,  Sy  =  1-01  and  a  =  10'^  S/m)  does  not  closely  follow  the  three  empirical 

models,  i.e.,  Tewari’s  model,  the  Jansky  and  Bailey  model  or  the  FITFl-R  model. 

Generally,  for  denser  foliage  (i.e.  Sr  >  1.05  and  a  >  10'^  S/m),  both  Tewari’s  model  and 

the  Jansky  and  Bailey  model  are  better  approximations  to  the  simulation  results  than  the 

FITFl-R  empirical  model  for  receiver  situated  at  higher  heights.  Generally,  the  FITFl-R 

model  gives  a  better  approximation  to  the  simulation  results  for  both  transmitter  and 
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receiver  situated  near  the  ground.  While  the  simulated  path  losses  do  not  exactly  match 
those  given  by  the  three  empirical  models,  the  difference  between  the  simulated  result 
and  empirical  result  is  approximately  10  dB  in  most  cases.  This  could  be  due  to  various 
factors.  While  the  empirical  models  were  proposed  to  be  an  approximation  to  the  loss 
experienced  for  foliage  propagation,  it  should  be  noted  that  these  models  were  formulated 
based  on  measured  data  from  foliage  in  a  particular  geographical  location.  Given  the 
lack  of  knowledge  on  the  electrical  properties  and  a)  of  the  actual  foliage,  assumed 
values  of  the  electrical  properties  were  used  in  the  simulations.  Hence,  it  was  expected 
that  there  would  be  difference  between  the  simulated  results  and  empirical  figures. 
Nevertheless,  given  the  observed  difference,  the  simulation  models  can  only  be 
considered  as  a  rough  approximation  to  actual  foliage  propagation. 

B,  COVERAGE  DIAGRAMS 

In  this  section,  the  effects  of  the  forest  on  the  coverage  of  a  transmitting  antenna 
immersed  in  foliage  are  examined.  The  coverage  diagram  for  a  transmitting  antenna  can 
be  obtained  from  the  electric  field  plots  of  a  CST  simulation.  A  simulation  is  first  carried 
out  for  the  transmitting  antenna  in  free  space.  This  is  to  establish  a  baseline  for 
comparison  purposes.  Subsequent  simulation  runs  are  carried  out  to  examine  the  effects 
of  the  forest  on  radiowave  propagation  with  respect  to  height  of  the  transmitter  in  foliage 
and  forest’s  electrical  properties.  The  general  observations  for  the  simulation  results  are 
discussed  at  the  end  of  this  section. 

1.  Transmitter  in  Free  Space 

A  dipole  antenna  is  set  up  as  the  transmitting  antenna.  The  antenna’s  height  is 
very  close  to  the  ground.  The  ground  is  represented  by  the  x-y  plane  ( =  0)  with  the 

boundary  condition  configured  as  ''electric  (Et  =  0)".  The  simulation  is  carried  out  with 
the  antenna  at  four  different  heights.  The  heights  ht  are  2  m,  5  m,  10  m  and  20  m.  The 
values  used  to  set  up  the  test  case  in  CST  are  summarized  in  Table  8.  Due  to  practical 
reasons,  an  adjacent  volume  of  space  is  defined  as  free  space  to  define  the  computational 
domain  for  the  simulation  of  fields  in  free  space.  This  is  achieved  via  the  background 
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properties  in  CST,  and  the  settings  are  shown  in  Figure  66.  E-field  and  Power  monitors 
are  added  to  monitor  the  fields  at  a  frequency  of  50  MHz. 


Table  8.  Values  used  to  set-up  dipole  antenna  transmitting  in  free  space. 


Parameter 

Description 

Value 

1 

n 

Length  of  dipole  antenna 

2.6  m 

2 

f 

Frequency  of  interest 

50  MHz 

3 

xd 

x-coordinate  of  dipole  center 

50  m 

4 

yd 

y-coordinate  of  dipole  center 

Om 

5 

zd 

z-coordinate  of  dipole  center 

2,  5,  10  or  20  m 

Figure  66.  Setting  for  background  properties  (for  coverage  diagram  test  cases). 


The  generated  E-fields  for  the  dipole  antenna  at  heights  of  2  m,  5  m,  10  m  and  20 
are  shown  in  Figures  67,  68,  69  and  70,  respectively. 
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Figure  67.  Coverage  diagram  for  a  dipole  antenna  in  free  space  (height  2  m). 

In  Figure  67,  the  transmitting  antenna  is  at  a  height  of  2  m.  It  can  be  observed 
that  there  is  only  one  main  lobe  radiating  out  from  the  antenna. 


Figure  68.  Coverage  diagram  for  a  dipole  antenna  in  free  space  (height  5  m). 


In  Figure  68,  the  transmitting  antenna  is  at  a  height  of  5  m.  It  can  be  observed 
that  there  are  two  main  lobes  radiating  out  from  the  antenna.  One  of  the  lobes  is 
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seemingly  along  the  ground,  while  the  other  is  at  an  angle  of  approximately  30°  from  the 
ground. 
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Figure  69.  Coverage  diagram  for  a  dipole  antenna  in  free  space  (height  10  m). 


In  Figure  69,  the  transmitting  antenna  is  at  a  height  of  10  m.  Compared  to  Figure 
28,  it  can  be  observed  that  the  number  of  lobes  increases  as  the  antenna  height  increases. 
In  this  case,  three  main  lobes  can  be  observed. 


Figure  70.  Coverage  diagram  for  a  dipole  antenna  in  free  space  (height  20  m). 
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In  Figure  70,  the  transmitting  antenna  is  at  a  height  of  20  m.  Compared  to  Figures 
28  and  29,  there  are  more  lobes  and  the  lobes  are  also  narrower. 

2,  Transmitter  in  Forest  (s^  =  1.065  and  a  =  0.000135  S/m) 

In  this  test  ease,  the  objeetive  is  to  examine  the  eoverage  diagram  with  the 
transmitting  dipole  antenna  immersed  in  foliage.  The  antenna  is  set  up  in  the  same 
manner  as  the  free  space  test  case.  A  rectangular  dielectric  block  is  added  to  the  model. 
The  antenna  and  forest  are  created  using  the  values  shown  in  Table  9.  Free  space  is 
defined  around  the  dielectric  block  using  the  settings  in  Figure  66.  The  amount  of  free 
space  added  is  50  m  above  the  forest  (i.e.,  along  the  positive  z-direction)  and  50  m  to  the 
right  of  the  forest  (i.e.,  along  the  positive  x-direction).  The  constructed  model  in  CST  is 
shown  in  Figure  7 1 . 


Table  9.  Values  used  to  set  up  model  for  examining  coverage  diagram  with 

transmitting  antenna  immersed  in  foliage. 


Parameter 

Description 

Value 

1 

n 

Length  of  dipole  antenna 

2.54  m 

2 

f 

Frequency  of  interest 

50  MHz 

3 

H 

Height  of  dielectric  block 

25  m 

4 

L 

Length  of  dielectric  block 

300  m 

5 

W 

Width  of  dielectric  block 

100  m 

6 

xd 

x-coordinate  of  dipole  center 

50  m 

7 

yd 

y-coordinate  of  dipole  center 

Om 

8 

zd 

z-coordinate  of  dipole  center 

2,  5,  10  or  20  m 

9 

6r 

Relative  permittivity 

1.065 

10 

a 

Conductivity 

0.000135  S/m 
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Figure  71.  Dipole  antenna  immersed  in  dielectric  forest  block  and  free  space  outside 

the  block. 


The  simulation  is  executed  for  varying  antenna  heights  of  2,  5  ,  10  and  20  m  with 
the  corresponding  E-field  plots  (generated  by  CST)  shown  in  Figures  72,  73,  74  and  75, 
respectively. 
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Figure  72.  Coverage  diagram  for  dipole  antenna  in  foliage  (height  2  m). 
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Comparing  Figure  72  and  Figure  67,  we  observe  that  the  addition  of  the  forest 
dielectric  block  causes  the  main  lobe  (i.e.  maxima)  along  the  ground  to  be  shorter.  For  a 
point  located  at  a  separation  distance  of  50  m  from  the  transmitter  and  just  above  the 
ground,  the  E-field  strength  is  —24.5  dBV/m  in  Figure  72  versus  —13.7  dBV/m  in  Figure 
67.  This  is  due  to  the  attenuation  as  the  EM  wave  propagates  through  the  lossy  dielectric 
media. 
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Eigure  73.  Coverage  diagram  for  dipole  antenna  in  foliage  (height  5  m). 

Comparing  Eigure  73  with  Eigure  68,  we  observe  that  there  are  still  two  main 
lobes  in  the  coverage  diagram.  However,  for  the  case  of  the  transmitter  immersed  in 
foliage,  the  maximum  of  each  lobe  is  reduced.  At  a  point  given  by  the  coordinates  x  = 
100  m,  y  =  0  m  and  z  =  35  m  just  above  the  foliage  block  along  one  of  the  maxima,  the 
E-field  strength  is  —25.1  dBV/m  in  Eigure  73  versus  —15.4  dBV/m  in  Eigure  68.  This  is 
because  the  EM  waves  experience  a  higher  attenuation  when  travelling  through  the  lossy 
dielectric  media  as  compared  to  free  space. 
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Figure  74.  Coverage  diagram  for  dipole  antenna  in  foliage  (height  10  m). 

Comparing  Figure  74  with  Figure  69,  we  observe  that  there  are  still  two  main 
lobes  in  the  eoverage  diagram.  However,  for  the  ease  of  the  transmitter  immersed  in 
foliage,  the  maximum  of  eaeh  lobe  is  reduced.  At  a  point  given  by  the  coordinates  x  = 
100  m,  y  =  0  m  and  z  =  38  m  above  the  foliage  block  along  one  of  the  maxima,  the  E- 
field  strength  is  —24.7  dBV/m  in  Figure  74  versus  —15.6  dBV/min  Figure  69.  This  is 
because  the  EM  waves  experience  a  higher  attenuation  when  travelling  through  the  lossy 
dielectric  media  as  compared  to  free  space. 

In  Figure  75,  the  main  lobes  can  be  faintly  observed.  While  the  number  of  main 
lobes  is  comparable  to  the  free  space  case  (Figure  70),  the  maxima  of  the  main  lobes  are 
greatly  reduced.  At  a  point  given  by  the  coordinates  x  =  100  m,  y  =  0m  and  z  =  28  m 
above  the  foliage  block  along  one  of  the  maxima,  the  E-field  strength  is  —22.4  dBV/m  in 
Eigure  75  versus  —15.4  dBV/m  in  Eigure  70.  The  lobes  are  due  to  the  reflections  of  EM 
waves  by  the  ground.  As  the  antenna  is  sited  higher,  the  path  length  that  the  reflected 
wave  has  to  travel  becomes  longer.  When  the  transmitter  is  sited  in  the  lossy  dielectric 
media,  the  reflected  waves  experience  higher  attenuation  as  compared  to  the  free  space 
case.  This  results  in  much  shorter  lobes  (as  evidenced  in  Eigure  75). 
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Figure  75.  Coverage  diagram  for /i/2  dipole  antenna  in  foliage  (height  20  m). 


3,  Transmitter  in  Forest  {Sr  =  1.15  and  a  =  0.00015  S/m) 

In  this  test  case,  the  coverage  diagram  is  examined  with  higher  e,-  and  a  values, 
which  represent  denser  foliage.  The  antenna  is  set  up  in  the  same  manner  as  the  free 
space  test  case.  A  rectangular  dielectric  block  is  added  to  the  model.  The  antenna  and 
forest  are  created  using  the  values  shown  in  Table  10.  Free  space  is  defined  around  the 
dielectric  block  using  the  settings  in  Figure  76.  The  amount  of  free  space  added  is  50  m 
above  the  forest  (i.e.,  along  the  positive  z-direction)  and  50  m  to  the  right  of  the  forest 
(i.e.,  along  the  positive  x-direction).  The  constructed  model  in  CST  is  similar  to  that 
shown  in  Figure  7 1 . 

The  simulation  is  executed  for  varying  antenna  heights  of  2,  5  ,  10  and  20  m  with 
the  corresponding  E-field  plots  (generated  by  CST)  shown  in  Figures  76,  77,  78  and  79, 
respectively. 


73 


Table  10.  Values  used  to  set  up  model  for  examining  eoverage  diagram  with 

transmitting  antenna  immersed  in  foliage. 


Parameter 

Description 

Value 

1 

£2 

Length  of  dipole  antenna 

2.54  m 

2 

f 

Frequency  of  interest 

50  MHz 

3 

H 

Height  of  dielectric  block 

25  m 

4 

L 

Length  of  dielectric  block 

300  m 

5 

W 

Width  of  dielectric  block 

100  m 

6 

xd 

x-coordinate  of  dipole  center 

50  m 

7 

yd 

y-coordinate  of  dipole  center 

Om 

8 

zd 

z-coordinate  of  dipole  center 

2,  5,  10  or  20  m 

9 

Relative  permittivity 

1.15 

10 

a 

Conductivity 

0.00015  S/m 

Comparing  Figure  76  and  Figure  72,  we  observe  that  the  main  lobe  (i.e.,  maxima) 


along  the  ground  is  slightly  shorter  in  the  denser  foliage  bloek.  For  a  point  loeated  at  a 
separation  distance  of  50  m  from  the  transmitter  and  just  above  the  ground,  the  E-field 
strength  is  —25  dBV/m  in  Figure  72  versus  —24.5  dBV/m  in  Figure  72.  This  is  due  to  the 
increased  attenuation  of  the  wave  in  the  denser  foliage  with  higher  values  of  Sr  and  a. 
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Figure  76.  Coverage  diagram  for  a  dipole  antenna  in  foliage  (height  2  m). 
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Figure  77.  Coverage  diagram  for  a  dipole  antenna  in  foliage  (height  5  m). 


Comparing  Figure  77  with  Figure  73,  we  observe  that  there  are  still  two  main 
lobes  in  the  coverage  diagram.  However,  the  maximum  of  each  lobe  is  reduced  in  the 
denser  foliage  block.  At  a  point  given  by  the  coordinates  x  =  100  m,  y  =  0m  and  z  =  35 
m  above  the  foliage  block  along  one  of  the  maxima,  the  E-field  strength  is  —25.9  dBV/m 
in  Figure  77  versus  —25.1  dBV/m  in  Figure  73.  This  is  due  to  the  increased  attenuation 
of  the  wave  in  the  denser  foliage  with  higher  values  of  Sr  and  a. 
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Figure  78.  Coverage  diagram  for  a  dipole  antenna  in  foliage  (height  10  m). 
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Comparing  Figure  78  with  Figure  74,  we  observe  that  there  are  still  two  main 
lobes  in  the  coverage  diagram.  However,  the  maximum  of  each  lobe  is  reduced  for  the 
case  of  the  denser  foliage  block.  At  a  point  given  by  the  coordinates  x=  100  m,  y  =  0m 
and  z  =  38  m  above  the  foliage  block,  the  E-field  strength  is  —27.3  dBV/m  in  Figure  78 
versus  —24.7  dBV/m  in  Figure  74.  This  is  due  to  the  increased  attenuation  of  the  wave  in 
the  denser  foliage  with  higher  values  of  e,-  and  a.  It  is  also  observed  that  the  point  (x  = 
100  m,  y  =  0  m,  z  =  38  m)  does  not  lie  exactly  on  the  maxima.  From  Figure  74,  the  E- 
field  strength  at  point  (x  =  101  m,  y  =  0  m,  z  =  33  m)  along  one  of  the  maxima  is  —25 
dBV/m. 
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Eigure  79.  Coverage  diagram  for  a  dipole  antenna  in  foliage  (height  20  m). 


In  Eigure  79,  the  difference  between  the  maxima  and  nulls  is  much  smaller. 
Hence,  the  lobes  in  Eigure  79  are  not  easily  discernible.  Comparing  Eigure  79  with 
Eigure  75,  we  observe  that  the  number  of  lobes  is  the  same.  However,  the  maximum  of 
each  lobe  is  reduced  for  the  case  of  the  denser  foliage  block.  At  a  point  given  by  the 
coordinates  x  =  100  m,  y  =  0m  and  z  =  28  m  above  the  foliage  block,  the  E-field  strength 
is  —28.3  dBV/m  in  Eigure  79  versus  —22.4  dBV/m  in  Eigure  75.  This  is  due  to  the 
increased  attenuation  of  the  wave  in  the  denser  foliage  with  higher  values  of  Sr  and  a. 
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4. 


Observations 


In  general,  it  is  observed  that  an  inerease  in  the  height  of  the  transmitter  inside  the 
foliage  results  in  more  maxima  (or  lobes).  In  addition,  the  angular  difference  between  the 
maxima  and  nulls  becomes  smaller  as  the  antenna  height  is  increased  until  near  the 
canopy  of  the  foliage.  As  the  transmitter  is  sited  higher  up  in  foliage,  the  path  length  of 
the  reflected  wave  in  foliage  becomes  longer  and  results  in  a  lower  E-field  strength.  For 
a  transmitter  sited  inside  a  denser  forest  with  higher  Sr  and  a  values,  the  lobes  are 
shortened  as  compared  to  ones  sited  in  a  less  dense  forest.  This  is  due  to  the  increased 
attenuation  in  denser  foliage,  which  is  a  result  of  higher  Sr  and  a  values. 

C.  COMMENTS  ON  SIMULATION  USING  CST 

General  comments  on  using  CST  for  the  simulation  of  radiowave  propagation  in 
foliage  are  made  in  this  section.  During  the  conduct  of  the  various  simulation  runs,  we 
made  several  observations  on  the  limitations  of  CST  with  respect  to  foliage  simulation. 
The  simulations  are  run  on  a  personal  computer  (PC)  equipped  with  Intel  Xeon  Quad- 
Core  2.53  GHz  processor  and  4  GB  RAM  memory.  The  PC  is  running  on  the  Windows  7 
Professional  operating  system.  One  main  limitation  on  the  simulation  with  respect  to 
total  simulation  time  is  the  maximum  frequency  used.  In  order  for  a  simulation  run  to  be 
completed  within  12  hours,  the  frequency  used  should  not  exceed  75  MHz  for  a  dielectric 
block  of  size  200  m  by  100  m  by  25  m.  An  increase  in  the  maximum  frequency  results  in 
an  exponential  increase  in  the  number  of  mesh  cells  generated  by  CST.  A  larger  number 
of  mesh  cells  requires  more  computer  resources  and  leads  to  an  increase  in  simulation 
time. 

D.  SUMMARY 

In  summary,  the  feasibility  of  using  a  single  dielectric  block  to  model  foliage 
propagation  at  the  HF/VHF  frequency  band  was  examined.  It  was  found  that  the 
proposed  simulation  model  using  CST  provides  a  rough  approximation  to  actual  foliage 
propagation.  From  the  simulation  results,  it  was  also  shown  that  the  path  loss  in  foliage 
is  affected  by  transmitter  height,  receiver  height  and  foliage  electrical  characteristic  {sr 
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and  a).  The  effect  of  foliage  on  a  transmitter  (located  inside  the  foliage)  was  also 
examined.  Compared  to  free  space,  the  lobes  from  a  transmitter  in  foliage  is  shorter  due 
to  the  increased  attenuation  experienced  by  radiowave  propagating  in  foliage.  Increasing 
the  height  of  a  transmitter  in  foliage  leads  to  a  smaller  difference  between  the  maxima 
and  nulls  in  the  coverage  diagram.  In  a  practical  application  such  as  radar,  this  would 
improve  consistency  in  tracking  a  target  flying  at  a  constant  height.  However,  the 
detection  range  is  also  reduced  due  to  the  increased  attenuation  through  foliage. 
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V.  SUMMARY  AND  CONCLUSION 


A,  SUMMARY 

The  main  objective  of  this  thesis  was  to  investigate  the  various  simulation  models 
for  radiowave  propagation  in  foliage  environments.  The  simulation  models  are  built 
upon  the  concept  of  representing  the  forest  as  a  single  dielectric  block.  The  model  is 
valid  for  frequencies  up  to  100  MHz.  There  are  two  important  parameters  that  represent 
the  forest  electrical  characteristics.  They  are  Sr  and  a.  These  parameters  affect  the 
radiowave  propagation  in  foliage.  The  values  used  for  Sr  and  a  in  the  simulation  models 
are  obtained  from  literature  that  document  the  values  derived  from  measured  data  in  real 
forest. 

The  EM  simulation  application,  CST  Studio  Suite,  is  used  for  the  simulation. 
After  the  simulation  is  completed,  MATLAB  was  used  to  process  the  field  data  and 
display  the  loss-distance  curves.  The  contour  plots  for  the  E-field  strength  are  given 
directly  by  CST. 

In  order  to  examine  the  suitability  of  the  proposed  simulation  models,  three 
empirical  models  are  selected  for  comparison  with  the  simulated  results.  The  models  are 
Tewari’s  model,  the  Jansky  and  Bailey  model  and  the  EITU-R  model.  The  effect  of  the 
foliage’s  characteristic  {Sr  and  a),  height  of  transmitter  and  height  of  receiver  were  also 
examined  in  this  thesis. 

The  coverage  diagram  for  a  transmitter  located  inside  foliage  was  also  examined 
in  this  thesis.  Different  values  for  the  foliage  electrical  characteristics  and  height  of 
transmitter  were  used  in  the  simulation  runs  to  examine  their  effects  on  the  coverage 
diagram 

B,  CONCLUSIONS 

The  results  and  data  analysis  from  the  simulation  show  that  the  proposed 
simulation  models  provide  a  rough  approximation  to  radiowave  propagation  in  an  actual 
rainforest  environment.  It  is  feasible  to  model  foliage  using  EM  simulation  applications. 
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From  the  simulation  technique  proposed  in  this  thesis,  an  RF  engineer  could  model  an  RF 
system  (especially  the  antenna)  that  is  being  developed  for  operations  in  foliage  in  CSX. 
The  performance  of  the  system  could  then  be  examined  during  the  design  phase.  This 
would  allow  the  designer  to  change  the  design  before  the  actual  system  was  built  and 
available  for  actual  field  tests. 

Based  on  the  simulation  results,  it  can  be  concluded  that  different  values  of  Sr,  a, 
height  of  transmitter  and  height  of  receiver  affect  the  path  loss  in  foliage.  The  coverage 
diagram  for  a  transmitter  immersed  in  foliage  is  also  similarly  affected  by  a  and  the 
height  of  the  transmitter.  For  a  transmitter  located  near  the  foliage  canopy,  the  resulting 
coverage  diagram  has  maxima  and  nulls  closer  to  each  other. 

C.  FUTURE  WORK 

1.  Model  the  Forest  as  Multiple  Smaller  Blocks 

Further  investigation  work  should  be  carried  out  to  determine  the  usefulness  of 
higher  fidelity  physical  model  of  the  forest.  Instead  of  representing  the  forest  as  a  single 
dielectric  block,  the  forest  could  be  decomposed  into  smaller  blocks.  This  “small  block” 
structure  would  be  useful  for  conducting  foliage  propagation  studies  at  a  higher 
frequency,  where  the  smaller  wavelength  is  considered  to  be  significant  compared  to  the 
dimensions  of  the  leaves  and  branches.  The  forest  could  also  be  further  sub-divided  into 
multiple  layers  to  represent  the  undergrowth  layer,  tree  trunks  and  canopy  layer. 
Different  Sr  and  a  could  be  assigned  to  the  different  layers. 

As  actual  trees  in  a  forest  are  not  uniform  in  height,  the  smaller  blocks  could  have 
different  heights.  The  heights  are  randomly  generated.  This  is  a  model  of  the  roughness 
of  the  air-canopy  interface.  Simulation  could  be  conducted  to  investigate  the  effect  of  the 
roughness  (i.e.,  non-uniform  canopy  heights)  on  the  propagation  of  lateral  waves  along 
the  air-canopy  interface. 


80 


2. 


Perform  the  Simulation  using  FEKO 


Another  possibility  for  future  work  is  to  eonduet  the  simulation  with  another  EM 
simulation  applieation  such  as  FEKO.  FEKO  is  an  EM  simulation  software  tool,  based 
on  computational  electromagnetics  (OEM)  techniques.  The  objectives  would  be  to 
determine  the  suitability  of  other  EM  simulation  application  for  conducting  propagation 
studies  in  foliage  environment  as  well  as  to  validate  the  results  obtained  in  CSX.  Another 
objective  is  to  also  find  an  application  which  could  simulate  foliage  propagation  at  higher 
frequencies  in  a  more  efficient  manner  than  CSX. 
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APPENDIX  A  MATLAB  CODE 


This  MATLAB  code  was  used  to  plot  the  loss-distanee  curves  for  the  empirieal 
models,  i.e.,  Tewari’s  model,  the  Jansky  and  Bailey  model  and  the  LITU-R  model. 


%Plot  Empirieal  Models  of  Propagation  in  Foliage 

clear  all 
close  all 

%Tewari's  model  model 
%Vertical  Polarization 
f  =  [50  200  500  800]; 

A  =[0  0.4989  0.3658  0.2661]; 
alpha  =  [0  0.0125  0.0135  0.0140]; 

B  =  [1.917  1.8358  0.9040  0.5331]; 
e  =  [V,  ’g’,  ’r’,  ’m']; 
set(goa,'YDir','reverse'); 

for  x=l:4 
for  i  =  10:200 

Lb(i-9)=-27.57+20*logl0(f(x))-20*logl0((A(x)*exp(-alpha(x)*i)/i)+(B(x)/(i^2))); 

d(i-9)=i; 

end 

plot(d,Lb,e(x)); 

set(gca,'YDir','reverse','FontName','Times  New  RomanVFontSize',  1 8); 
hold  all; 

xlabel('Distanoe,  d  (m)VFontName','Times  New  Roman','FontSize',18) 
ylabel('Loss,  L  b  (dB)','FontName','Times  New  Roman','FontSize',18) 
end 

legend(’50  MHzV200  MHzV500  MHz',’800  MHz'); 

clear  d; 
clear  f; 

%Jansky  and  Bailey  model 
f=[25  50  100  250  400]; 

A  =[0  0  0.615  0.759  1.02]; 
alpha  =[0  0  0.045  0.050  0.055]; 

B  =  [0.00212  0.00106  0.000529  0.000443  0.000523]; 

c  =  ['b',  'g',  'r',  ’m’,’k’]; 

set(gca,'YDir','reverse'); 

figure; 
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for  x=l:5 
for  i  =  10:200 
z  =  i/1609; 

LJ(i-9)=36.57+20*logl0(f(x))-20*logl0((A(x)*exp(- 
1 609 *  alpha(x)  *  z)/z)+(B(x)/(z^2))) ; 
d(i-9)=i; 
end 

plot(d,LJ,c(x)); 

set(gca,'YDirVreverseVFontNameVTimes  New  RomanVFontSize',18); 
hold  all; 

xlabel(’Distanoe,  d  (m)VFontNameVTimes  New  RomanVFontSize',18) 
ylabel('Loss,  L  b  (dB)VFontNameVTimes  New  RomanVFontSize',18) 
end 

legend(’25  MHz’/SO  MHz'/lOO  MHz', '250  MHz', '400  MHz'); 

elear  d; 
elear  f; 

%Lateral  ITU-R  (LITU-R)  model 
%VHF  bands 
elear  d; 
elear  f; 

f=[30  50  100  300]; 

T  =  5; 

R  =  5; 

e  =  ['b'  'g'  'r'  'm']; 
figure; 
for  X  =  1:4 
for  i  =  20:200 

Litur(i-19)=40*logl0(i)-20*logl0(T)-20*logl0(R)+0.48*f(x)^0.43*H0.13; 

d(i-19)=i; 

end 

plot(d,Litur,c(x)); 

set(gea,'YDir','reverse','FontName','Times  New  Roman','FontSize',18); 
hold  all; 
end 

xlabel('Distanee,  d  (m)','FontName', 'Times  New  Roman','FontSize',18) 
ylabel('Loss,  L  (dB)','FontName','Times  New  Roman','FontSize',18) 
legend('30  MHz  H_t=H_r=5m','50  MHz  H_t=H_r=5m','100  MHz  H_t=H_r=5m','300 
MHz  H_t=H_r=5m'); 

xlim([0  200]); 
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This  MATLAB  code  was  used  to  plot  the  loss-distance  curves  for  the  simulated 
results  in  Chapter  IV. 


%  Plot  Loss  versus  Seperation  Distanee  Curves 

%  Using  ASCII  fde  of  powerflow  data  (absolute  values)  exported  from  CSX  Studio 
%  Only  for  fields  projeeted  onto  the  X-Z  plane  in  CSX 
%  ereated  by  Chan  Chung  Wei,  ECE  Dept,  Naval  Postgradute  Sehool 

clear  all 
elose  all 


%  Import  data  from  Aseii  file  into  Matlab 
filename  =  input('Enter  file  name: 
file  =  importdata(filename); 
len  =  length(fde.data); 

L  =  input('max  x  eoord  of  block  (len):  '); 

W  =  input('max  y  coord  of  block  (wt):  '); 

H  =  input('max  z  coord  of  block  (ht):  ’); 

EO  =  input('x  origin  of  block:  '); 

WO  =  input('y  origin  of  block:  '); 

HO  =  input('z  origin  of  block:  '); 

xO  =  input('dipole  centre,x:  ’); 

yO  =  input('dipole  centre,y:  '); 

zO  =  input('dipole  centre,z:  '); 

zoffset  =  input('Enter  offset  of  z  steps  in  data  file:  ’); 

%offet  obtained  by  looking  at  z  data  in  ASCII  file 

%  e.g.  offset  for  z  =  3.1  is  0.1 


%  Caleulate  separation  distanee  and  poynting  veetor  at  eaeh  point  (x,y,z) 

j  =  0; 

for  i=l:len-E0 

x(i)=single(file.data(i,l));  %import  x 
y(i)=single(file.data(i,2));  %import  y 
z(i)=single(fde.data(i,3));  %import  z 
PxRe(i)=single(fde.data(i,4));  %import  Px  Re 

%import  Py  Re 
%import  Pz  Re 
%import  Px  Imag 
%import  Py  Imag 
%import  Pz  Imag 


PyRe(i)=single(fde.data(i,5)); 

PzRe(i)=single(file.data(i,6)); 

Pxlm(i)=single(fde.data(i,7)); 

Pylm(i)=single(fde.data(i,8)); 

Pzlm(i)=single(file.data(i,9)); 


if  (x(i)<=L  &&  x(i)>=x0)  &&  (y(i)<=W  &&  y(i)>=W0)  &&  (z(i)<=H  &&  z(i)>=H0) 

j=j+i; 
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r(j)  =  sqrt((x(i)-x0)^2+(y(i)-y0)^2+(z(i)-z0)^2);  %calculate  distance 
Px(j)  =  sqrt(PxRe(i)^2+PxIm(i)^2); 

Py(j)  =  sqrt(PyRe(i)^2+PyIm(i)^2); 

Pz(j)  =  sqrt(PzRe(i)^2+PzIm(i)^2); 

Pr(j)  =  -10*logl0(sqrt(Px(j)^2+Py(j)^2+Pz(j)^2));  %calculate  magnitude  of 
Poinyting  vector 

%loss(j)  =  0-Pr(j);  Transmit  Power  is  1  W  =  0  dBW 
htG)=z(i); 
end 
end 


hl=input(’lower  height  limit  to  plot: '); 
hu=input('upper  height  limit  to  plot:  ’); 

for  disph=hl:hu 
clear  dispr; 
clear  dispPr; 
clear  index; 

dht  =  disph+zoffset; 

j  =  0; 

for  i=l:length(ht) 
if  ht(i)==dht 

j=j+i; 

indexG)  =  i; 
end 
end 


for  i=l:j 

dispr(i)=r(index(i)); 

dispPr(i)=Pr(index(i)); 

end 

semilogx(dispr,dispPr); 

set(gca,'YDirVreverseVFontNameVTimes  New  Roman'/FontSize',  1 8); 
hold  all 
end 


title('Plot  of  Loss  vs  Distance  at  different  Height,  Hf ); 
xlabel('Distanee,  r  (m)') 
ylabel('Loss,  L  (dB)') 

%plot  Tewari  Basic  Transmission  Loss  Model 
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f=50; 

A  =  0; 
alpha  =  0; 

B=  1.917; 
for  1=  10:300 

Lb(i-9)=-27.57+20*logl0(f)-20*logl0((A*exp(-alpha*i)/i)+(B/(i^2))); 

d(i-9)=i; 

end 

semilogx(d,Lb,'— r'); 

%  Jansky  and  Bailey  Model 

clear  d; 

f=50; 

A  =  0; 
alpha  =  0; 

B  =  0.00106; 

for  i=  10:300 
z  =  i/1609; 

LJ(i-9)=36.57+20*logl0(f)-20*logl0((A*exp(-1609*alpha*z)/z)+(B/(z^2))); 

d(i-9)=i; 

end 

semilogx(d,LJ,'— m'); 

%lateral  ITU-R  model 

clear  d; 

f=50; 

%  plot  curve  for  receiver  height  from  hu  to  hi 
for  X  =  l:hu-hl+l 
hr  =  hl-l+x+offset; 
for  i  =  20:300 

Litur(i-19)=40*logl0(i)-20*logl0(z0)-20*logl0(hr)+0.48*R0.43*i^0.13; 

d(i-19)=i; 

end 

semilogx(d,Litur,'- . ') 
end 

xlim([20  150]); 

%legend('Ht  =  20VHt  =  21', 'Ht  =  22', 'Ht  =  23','Ht  =  24',  'Ht  =  25','Tewari  model', 'Jansky 
&  Bailey  Model','Lateral  ITU-R'); 
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